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O' ABSTRACT 
I We present the PPAK Integral Field Spectroscopy (IFS) Nearby Galaxies Survey: 

\^ • PINGS, a 2-dimensional spectroscopic mosaicking of 17 nearby disk galaxies in the 

(-H ' optical wavelength range. This project represents the first attempt to obtain continu- 

Q^, ous coverage spectra of the whole surface of a galaxy in the nearby universe. The final 

I ' data set comprises more than 50000 individual spectra, covering in total an observed 

O I area of nearly 80 arcmin^. The observations will be supplemented with broad band 

and narrow band imaging for those objects without public available images in order 
^ ' to maximise the scientific and archival value of the data set. In this paper we describe 

the main astrophysical issues to be addressed by the PINGS project, we present the 
galaxy sample and explain the observing strategy, the data reduction process and all 
uncertainties involved. Additionally, we give some scientific highlights extracted from 
^ I the first analysis of the PINGS sample. A companion paper will report on the first 

■ results obtained for NGC 628: the largest IFS survey on a single galaxy. 

0^ ' Key words: Surveys - methods: observational - techniques: spectroscopic - galaxies: 

general - galaxies: abundances - ISM: abundances 
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1 INTRODUCTION 2D spectroscopic mosaics of a representative sample of 

nearby spiral galaxies, using the unique instrumental ca- 
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by means of pixel-resolved maps across the disks of nearby 
galaxies. These spectral maps will allow us to test, confirm, 
and extend the previous body of results from small-sample 
studies, while at the same time open up a new frontier of 
studying the two-dimensional metallicity structure of disks 
and the intrinsic dispersion in metallicity. Furthermore, the 
large body of data arising from these studies will also allow 
us to test and strengthen the diagnostic methods that are 
used to measure H II region abundances in galaxies. 

Previous works have used multi-object instruments to 
obtain simultaneous spectra of H II regions in a disk galaxy 
(e.g. Roy & Walsh 1988; Kennicutt & Garnett 1996; Mous- 
takas & Kennicutt 2006a), or narrow-band imaging of spe- 
cific fields to obtain information of star forming regions and 
the ionized gas (e.g. Scowen et al. 1996). One important at- 
tempt is represented by the SAURON project (Bacon et al. 
2001), which is based on a panoramic lenslet array spec- 
trograph with a relatively large FOV of 33x41 arcsec'^. 
SAURON was specifically designed to study the kinemat- 
ics and stellar populations of a sample of nearby elliptical 
and lenticular galaxies. The application of SAURON to spi- 
ral galaxies was restricted to the study of spiral bulges. A 
recent effort by Rosolowsky & Simon (2008) plans to obtain 
spectroscopy for ~ 1000 H II regions through the M 33 Metal- 
licity Project, using multi-slit observations. On the other 
hand, Blanc et al. (2009) obtained IPS observations of the 
central region of M51 (~ 1.7 arcmin^), using the VIRUS-P 
instrument. However, in spite of the obvious advantages of 
the IPS technique in tackling known scientific problems and 
in opening up new lines of research, 2D spectroscopy is a 
method that has been used relatively infrequently to study 
large angular-size nearby objects. Recently, PPAK was used 
successfully to map the Orion nebula, obtaining the chem- 
ical composition through strong line ratios (Sanchez et al. 
2007b). Likewise, PMAS in the lens-array configuration was 
used to map the spatial distribution of the physical proper- 
ties of the dwarf HII galaxy IIZw70 (Kehrig et al. 2008), 
although covering just a small FOV (~ 32 arcsec). 

Despite these previous efforts toward IPS of nearby 
galaxies, the application to obtain complete 2D information 
in galaxies is a novel technique. Reasons for the lack of stud- 
ies in this area include small wavelength coverage, fibre-optic 
calibration problems, but mainly the limited POV of the in- 
struments available worldwide. Most of these IPUs have a 
POV of the order of arcsec, preventing a good coverage of 
the target galaxies on the sky in a reasonable time, even 
with a mosaicking technique. Furthermore, in some cases the 
spectral coverage is not appropriate to measure important 
diagnostic emission-lines used in chemical abundance stud- 
ies. Moreover, the complex data reduction and visualisation 
imposes a further obstacle to more ambitious projects based 
on 2D spectroscopy. To our knowledge there has been no 
attempt to obtain point-by-point spectra over a large wave- 
length range of the whole surface of a galaxy covering all 
HII regions within it. Similarly to SAURON for early- type 
galax;ies, PINGS will provide the most detailed knowledge of 
star formation and gas chemistry across a late-type galaxy. 
This information is also relevant for interpreting the inte- 
grated colours and spectra of high redshift sources. In that 
respect, PINGS represents a leap in the study of the chem- 
ical abundances and the global properties of galaxies. 

The objectives of this paper are: 1) to provide the back- 



ground information of the PINGS survey; including a de- 
tailed description of the observations, all the data reduction 
techniques implemented (some of them novel in the treat- 
ment of IPS data), and all the uncertainties involved in this 
process; 2) to offer a general description of the procedures 
involved in IPS observations of this kind; and 3) to present 
the PINGS data products and the archival value of this sur- 
vey. The paper is organised as follows. In § 2 we describe 
the core scientific objectives of the PINGS project and we 
discuss some of the many applications of the data set. In 
§ 3 we describe the properties of the galaxy sample selected, 
while in § 4 we explain the observational strategy and the 
PINGS observations themselves. In § 5 we present the des- 
cription of the complex data reduction involved in this IFS 
survey, while in § 6 we describe the sources and magnitudes 
of the errors and uncertainties in the data sample. In § 7, we 
summarise the basic properties of the data, showing a few 
science case examples extracted from the data set, including 
the integrated properties, emission line maps and a compar- 
ison of line intensity ratios with previously published studies 
for some galaxies. Finally, in § 8 we give a summary of the 
article. 



2 SCIENTIFIC OBJECTIVES 

The study of chemical abundances in galaxies has been sig- 
nificantly benefited from the vast amount of data collected 
in recent years, either at the neighborhood of the Sun, or 
at high redshifts, especially on large scale surveys such as 
the SDSS or the 2dP galaxy redshift survey (CoUess et al. 
2001). Most studies derived from these observations have 
focused on linking the properties of high redshift galaxies 
with nearby objects, as an attempt to understand the prin- 
ciples of the formation and chemical evolution of the galaxies 
(e.g. see Pettini 2006). Historically, the metal content of low 
redshift galaxies has been determined through the nebular 
emission of individual HII regions at discrete spatial po- 
sitions, these measurements provide hints on the chemical 
evolution, stellar nucleosynthesis and star formation histo- 
ries of spiral galaxies. The chemical evolution is dictated by 
a complex array of parameters, including the local initial 
composition, the distribution of molecular and neutral gas, 
star formation history (SFH), gas infall and outflows, radial 
transport and mixing of gas within disks, stellar yields, and 
the initial mass function (IMF). Although it is difficult to 
disentangle the effects of the various contributors, measure- 
ments of current elemental abundances constrain the possi- 
ble evolutionary histories of the existing stars and galaxies, 
and thus the importance of the accurate determination of 
the chemical composition among different galaxy types. 

Different studies have shown a complex link between the 
chemical abundances of galaxies and their physical proper- 
ties. Such studies are only able to accurately measure the 
first two moments of the abundance distribution -the mean 
metal abundances of disks and their radial gradients- and on 
characterising the relations between these abundance prop- 
erties and the physical properties of the parent galaxies, for 
example galactic luminosity, stellar and dynamical mass, cir- 
cular velocity, surface brightness, colors, mass-to-light ra- 
tios, Hubble type, gas fraction of the disk, etc. These stud- 
ies have revealed a number of important scaling laws and 
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systematic patterns including luminosity-metallicity, mass- 
metallicity, and surface brightness vs. metallicity relations 
(e.g. Skillman et al. 1989; Vila-Costas & Edmunds 1992; 
Zaritsky et al. 1994; Tremonti et al. 2004), eflective yield 
vs. luminosity and circular velocity relations (e.g. Garnett 
2002), abundance gradients and the effective radius of disks 
(e.g. Diaz 1989), and systematic differences in the gas- phase 
abundance gradients between normal and barred spirals (e.g. 
Zaritsky et al. 1994; Martin & Roy 1994). However, these 
studies have been limited by the rmmber of objects sam- 
pled, the number of H II regions observed and the coverage 
of these regions within the galaxy surface. 

In order to obtain a deeper insight of the mechanisms 
that rule the chemical evolution of galaxies, we require the 
combination of high quality multi- wavelength data and wide 
field optical spectroscopy in order to increase significantly 
the number of HII regions sampled in any given galaxy. 
The PINGS project was conceived to tackle the problem 
of the 2-dimensional coverage of the whole galaxy surface. 
The imaging spectroscopy technique applied in PINGS pro- 
vides a powerful tool for studying the distribution of physical 
properties in nearby well-resolved galaxies. PINGS was spe- 
cially designed to obtain complete maps of the emission-line 
abundances, stellar populations, and reddening using an IFS 
mosaicking imaging, which takes advantage of what is cur- 
rently one of the world's widest FOV IFU. With this novel 
spectroscopic technique, the data can be used to derive: 1) 
oxygen abundance distributions based on a suite of strong- 
line diagnostics incorporating absorption-corrected Ha, H/3, 
[Oil], [OIII], [Nil], and [SII] line ratios; 2) local nebular red- 
dening estimates based on the Balmer decrement; 3) mea- 
surements of ionization structure in HII regions and dif- 
fuse ionized gas using the well-known and most updated 
forbidden-line diagnostics in the oxygen and nitrogen lines; 
4) rough fits to the stellar age mix from the stellar spectra. 

The resulting spectral maps and ancillary data will be 
used to address a number of important astrophysical issues 
regarding both the gas-phase and the stellar populations in 
galaxies. For example, one application will be able to test 
whether the metal abundance distributions in disks are ax- 
isymmetric. This is usually taken for granted in chemical 
evolution models, but one might expect strong deviations 
from symmetry in strongly lopsided, interacting, or barred 
galaxies, which are subject to large scale gas flows. Another 
important goal is to place strong limits on the dispersion 
in metal abundance locally in disks; there is evidence for a 
large dispersion in some objects such as NGC 925 or M 33 
(Rosolowsky & Simon 2008), but it is not clear from those 
data whether the dispersion is due to non-axisymmetric 
abundance variations, systematic errors in the abundance 
measurements, or a real local dispersion. Yet another by 
product of our analysis will be point by point reddening 
maps of the galaxies, which can be combined with UV, 
Hq, and infrared maps to derive robust, extinction-corrected 
maps of the SFR. 

PINGS can also provide a very detailed knowledge of 
the role played by star formation in the cosmic life of galax- 
ies. All the important scaling laws previously mentioned tell 
us that, once born, stars change the ionization state, the 
kinematics and chemistry of the interstellar medium and, 
thus, change the initial conditions of the next episode of star 
formation. Substantially, star formation is a loop mechanism 



which drives the luminosity, mass and chemical evolution of 
each galaxy (leaving aside external agents like interactions 
and mergers). The details of such a complex mechanism are 
still not well established observationally and not well devel- 
oped theoretically, and limit our understanding of galaxy 
evolution from the early universe to present day. In com- 
bination with ancillary data, the flux maps computed from 
the PINGS data will be used to study both the most recent 
star formation activity of the targets and the older stellar 
populations. We will be able to identify the gas and stellar 
features responsible for the observed spectra, to derive the 
dependence of the local star formation rate on the local sur- 
face brightness, a key recipe for modelling galaxy evolution 
and the environmental dependence of star formation. These 
data will also provide an important check for interpreting 
the integrated broad-band colours and spectra of high red- 
shift sources. 



3 SAMPLE DESCRIPTION 

In order to achieve the scientific goals described above, we 
incorporated a diverse population of galaxies, adopting a 
physically based approach to defining the PINGS sample. 
We decided to observe a set of local spiral galaxies which 
were representative of different galaxy types. However, the 
size and precise nature of the sample were heavily influ- 
enced by a set of technical considerations, the principal lim- 
iting factor being the FOV of the PPAK unit. We wanted to 
observe relatively nearby galaxies to maximise the physical 
linear resolution using the mosaicking technique. However, 
we also had to take into account the limitations imposed by 
the amount of non-secure observing time and meteorological 
conditions for the granted runs. Therefore, the sample size 
was dictated by a balance between achieving a representa- 
tive range of galaxies properties and practical limitations in 
observing time. 

When constructing the sample, we also took into ac- 
count a range of other galaxy properties, such as inclination 
(with preference to face-on galaxies), surface brightness, bar 
structure, spiral arm structure, and environment (i.e. iso- 
lated, interacting and clustered). We favoured galaxies with 
high surface brightness and active star formation so that 
we could have a good distribution of H II regions across the 
galaxies surface. Galaxies with bars and/or non- typical spi- 
ral morphology were also preferred. The flnal selection of 
galaxies also took into account practical factors such as op- 
timal equatorial right ascension and declination for the lo- 
cation of Calar Alto observatory, and the observable time 
per night for a given object above a certain airmass (due to 
problems of differential atmospheric refraction). 

The PINGS sample consists of 17 galaxies within a ma- 
ximum distance of 100 Mpc; the average distance of the sam- 
ple is 28 Mpc (for Ho = 73 kms~^ Mpc"^). The final sam- 
ple includes normal, lopsided, interacting and barred spirals 
with a good range of galactic properties and SF environ- 
ments with multi-wavelength public data. A good fraction 
of the sample belongs to the Spitzer Infrared Nearby Galax- 
ies Survey (SINGS, Kennicutt et al. 2003), which ensures a 
rich set of ancillary observations in the UV, infrared, H I and 
radio. 

The sample objects were given a different observing pri- 
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Figure 1. Digital Sky Survey images for a selection of objects included in the PINGS sample. On each image, an array of hexagonal 
fields corresponding to the PPAK aperture is superimposed on scale, showing the IFU mosaicking technique and the observed positions. 
All images are 10' X 10' and displayed in top-north, left-east configuration. 



ority based on the angular size of the objects, the num- 
ber of PPAK adjacent pointings necessary to complete 
the mosaic, and the scientific relevance of the galaxy. Ta- 
ble 1 gives a complete listing of the PINGS sample with 
some of their relevant properties. The first priority was as- 
signed to medium-size targets such as NGC 1058, NGC 1637, 
NGC 3310, NGC 4625 and NGC 5474 which are bright, face- 
on spirals of very different morphological type, with many 
sources of ancillary data and could be covered with relatively 
few IFU pointings. The second priority was given to smaller 
galaxies which fit perfectly in terms of size and acquisition 
time for the periods during the night when the first priority 
objects were not observable (due to a high airmass or bad 
weather conditions) and/or in the case their mosaicking was 
completed. 



NGC 628 (Messier 74) is a special object among the se- 
lected galaxies and the most interesting one of the sample. 
NGC 628 is a close, bright, grand-design spiral galaxy which 
has been extensively studied. With a projected optical size 
of 10.5 X 9.5 arcmin, it is the most extended object of the 
sample. Although it could be considered too large to be fully 
observed in a realistic time, we attempted the observation of 
this galaxy considering that the spectroscopic mosaicking of 
NGC 628 represents the real 2-dimensional scientific spirit 
of the PINGS project (see Sanchez et al. 2009, hereafter Pa- 
per II). Such a large galaxy would offer us the possibility to 
assess the body of results from the rest of the small galaxies 
in the sample and would allow us to study the 2D metal- 
licity structure of the disk and second order properties of 
the abundance distribution. The observations of this galaxy 
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spanned a period of three-years. Hitherto, NGC 628 repre- 
sents the largest area ever covered by an IFU, as described 
briefly in the next section and in detail in Paper II. A special 
priority was also given to NGC 3184, galaxy which falls be- 
tween the medium first-priority and large size galaxies. The 
observations for this object spanned for 2 years, obtaining 
an almost complete mosaicking (see Figure I), making this 
galaxy the 2nd largest object of the sample. 

Figure I shows Digital Sky Survey^ images for a se- 
lection of galaxies listed in Table 1. The mosaicking of the 
largest objects in the sample, NGC 628 and NGC 3184, con- 
sist of 34 and 16 individual IFU pointings respectively, cov- 
ering almost completely the spiral arms of these two bright 
grand-design galaxies. The outlying pointings of NGC 1058 
and the mosaicking configurations of NGC 3310, NGC 6643 
and NGC 7771 are explained in the next section. 

In summary, the PINGS sample was selected in a careful 
way to find a trade-off between the size of the galaxies, their 
morphological types, their physical properties and the prac- 
tical limitations imposed by the instrument and the amount 
of observing time. The result is a comprehensive sample of 
galaxies with a good range of galactic properties and avail- 
able multi-wavelength ancillary data, in order to maximise 
both the original science goals of the project and the possible 
archival value of the survey. 



4 OBSERVATIONS 

Observations for the PINGS galaxies were carried out at 
the 3.5m telescope of the Calar Alto observatory with the 
Postdam Multi Aperture Spectrograph, PMAS (Roth et al. 

2005) in the PPAK mode (Verheijen et al. 2004; Kelz et al. 

2006) , i.e. "a retrofitted bare fibre bundle IFU which ex- 
pands the FOV of PMAS to a hexagonal area with a foot- 
print of 74 X 65 arcsec, with a filling factor of 65% due to 
gaps in between the fibres" . The PPAK unit features a cen- 
tral hexagonal bundle with 331 densely packed optical fibres 
to sample an astronomical object at 2.7 arcsec per fibre. 
The sky background is sampled by 36 additional fibres dis- 
tributed in 6 mini-IFU bundles of 6 fibres each, in a circular 
distribution at ~ 90 arcsec of the centre and at the edges of 
the central hexagon. The sky-fibres are distributed among 
the science ones in the pseudo-slit, in order to have a good 
characterisation of the sky. Additionally, 15 fibres can be 
illuminated directly by internal lamps to calibrate the in- 
strument. 

All sample galaxies were observed using the same tele- 
scope and instrument set-up. We used the V300 grating, 
covering a wavelength between 3700-7100 A with a reso- 
lution of ~ 10 A FWHM, corresponding to ~ 600 kms~^. 
With this set-up we cover all the optical strong emission 
lines used in typical abundance diagnostic methods (Sanchez 
et al. 2007b). For the particular setup used in the PINGS 

^ The Digitized Sky Survey was produced at the Space Telescope 
Science Institute under U.S. Government grant NAG W-2166. 
The images of these surveys are based on photographic data ob- 
tained using the Oschin Schmidt Telescope on Palomar Mountain 
and the UK Schmidt Telescope. The plates were processed into 
the present compressed digital form with the permission of these 
institutions. 



survey there was no need to use a order separating filter. The 
main reason being that the efficiency of the instrument -|- 
telescope system (mostly the grating refiectivity and the fi- 
bre transmission), drops dramatically at wavelengths shorter 
than ~ 3600 A, where the transmission is 1/40000 of the one 
at the peak intensity for the V300 grating (~ 5400 A) and 
1/10000 of the value at redder wavelengths covered by our 
survey (~ 6500 A). At bluer wavelengths the efficiency is 
even lower. Therefore, the system by itself blocks any pos- 
sible 2nd order contamination up to ~ 7200 A, and only at 
larger wavelengths it is required an order separation filter 
(Kelz et al. 2006). The exposure times were calculated from 
previous experience with the instrument in order to obtain 
spectroscopy with S/N J5 20 in the continuum and S/N 
50 in the Ha emission line for the brightest H II regions with 
the given grating. 

Different observing strategies were implemented de- 
pending on the size and priority of the targets. For those ob- 
jects with relatively small angular size, single PPAK point- 
ings would not sample the surface of the galaxy with enough 
spatial resolution, due to the incomplete filling factor of the 
fibre bundle. In this dithering procedure was applied 

(Sanchez et al. 2007a). For each individual position in dither- 
ing mode, the first exposure was recorded and then, two 
consecutive exposures with the same acquisition time were 
recorded, but with small offsets of A(RA, Dec) = (1.56, 0.78) 
and (1.56,-0.78) arcsec with respect to the first exposure. 
The advantage of this method is that all gaps of the original 
exposure are covered, and every single point of the dithered 
field is spectroscopically sampled within the resolution. The 
pitfalls are that the exposure time and the amount of data 
to be processed is triple a normal frame, preventing the pos- 
sibility of applying this method to large mosaics. We used a 
mean acquisition time per PPAK field in dithering mode (in- 
cluding set-up + integration time) of 2 x 600 sec. per dither- 
ing position for a total of 60 min. exposure; and 3 x 600 sec. 
for non-dithered frames. 

The observations extended over a period of three-years 
with a total of 19 observing nights distributed during dif- 
ferent observing runs and seasons. For all the objects in 
the sample, the first exposure was centered in a given ge- 
ometrical position which, depending on the morphology or 
a particular mosaicking pattern, may or may not coincide 
with the bright bulge of the galaxy. Consecutive pointings 
followed in general a hexagonal pattern, adjusting the mo- 
saic pointings to the shape of the PPAK science bundle as 
shown in Figure 1. Due to the shape of the PPAK bundle 
and by construction of the mosaics, 11 spectra of each point- 
ing corresponding to one edge of the hexagon, overlap with 
the same number of spectra from the previous pointing. This 
pattern was selected to maximise the covered area, but to 
allow enough overlapping to match the different exposures 
taken under variable atmospheric conditions. Exceptions are 
NGC 2976, NGC 3310, NGC 6643 and NGC 7770 in which 
the mosaics were constructed to optimise the galaxy surface 
coverage as explained below. 

4.1 NGC 628 

NGC 628 (or M 74) is an extensively studied isolated grand- 
design Sc spiral galaxy at a distance of 9.3 Mpc in the con- 
stellation of Pisces. The observations for this galaxy totaled 
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Table 1. Galaxy properties of the PINGS sample. Col. (1): Galaxy name. Col. (2): Morphological type from the R3C 
catalog{de Vaucouleurs et al. 1991). Col. (3): Distances in Mpc, references: NGC628, Hendry et al. (2005); NGC1058, 
Eastman et al. (1996); NGC 1637, Saha et al. (2006); NGC 2976, Karachentsev et al. (2002); NGC 3184, Leonard et al. (2002); 
NGC 3310, Terry et al. (2002); NGC 5474, Drozdovsky & Karachentsev (2000); NGC 6643, Willick et al. (1997); NGC 4625, 
NGC 6701, NGC 7771, & Stephan's Quintet: Galactocentric GSR distances derived from the redshift, assumming a value of 
Ho = 73 kms-l Mpc-^. Col. (4): Projected size, major and minor axes at the B25 mag arcscc from R3C, except NGC 7318 
from Jarrett et al. (2003). Col. (5): Absolute B-band magnitude calculated from the apparent magnitude listed in the R3C 
catalog and the adopted distances to the system. Col. (6): Redshift, references: NGC 628, Lu et al. (1993); NGC 4625, 
Fisher et al. (1995); NGC 6701, Theureau et al. (1998); NGC 7770, Woods et al. (2006); NGC 1637, NGC 3310, Haynes 
et al. (1998); NGC 2976, NGC 5474, Falco et al. (1999, The Updated Zwicky Catalog); NGC 1058, NGC 3184, NGC 6643, 
NGC 7771, Springob et al. (2005); NGC 7317, NGC 7318a, NGC 7318b, Hickson et al. (1992); NGC 7319, NGC 7320, Nishiura 
et al. (2000). Col. (7): Heliocentric velocities calculated from v = zc, with no further correction applied. Col. (8): Galaxy 
inclination angle based on the B25 mag arcsec"^ from R3C. Col. (9): Galaxy position angle, measured positive NE, in the 
B25 mag arcsec-2 except for NGC 1058, NGC 1637, NGC 7317, NGC 7319, NGC 7770, which are based on the Ks-band 
(Jarrett et al. 2003), and NGC 3310, NGC 4625, NGC 5474 based on the r-SDSS band. Col. (10): Galaxy location. 



Object 


Type 


Distance 

(Mpc) 


Projected size 

(arcmin) 


Mb 


2 


V0 

(kms~^) 


i 


P.A. 


Constellation 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


NGC 628 


SA(s)c 


9.3 


10.5 X 9.5 


-19.9 


0.00219 


657 


24 


25 


Pisces 


NGC 1058 


SA(rs)c 


10.6 


3.0 X 2.8 


-18.3 


0.00173 


519 


21 


95 


Perseus 


NGC 1637 


SAB(rs)c 


12.0 


4.0 X 3.2 


-18.9 


0.00239 


717 


36 


33 


Eridanus 


NGC 2976 


SAc pec 


3.6 


5.9 X 2.7 


-16.9 


0.00008 


24 


63 


143 


Ursa Major 


NGC 3184 


SAB(rs)cd 


11.1 


7.4 X 6.9 


-19.9 


0.00194 


582 


21 


135 


Ursa Major 


NGC 3310 


SAB(r)bc 


17.5 


3.1 X 2.4 


-20.1 


0.00331 


993 


39 


163 


Ursa Major 


NGC 4625 


SAB(rs)m 


9.0 


2.2 X 1.9 


-16.9 


0.00203 


609 


29 


30 


C. Venatici 


NGC 5474 


SA(s)cd 


6.8 


4.8 X 4.3 


-17.9 


0.00098 


294 


27 


91 


Ursa Major 


NGC 6643 


SA(rs)c 


20.1 


3.8 X 1.9 


-19.8 


0.00495 


1485 


60 


37 


Draco 


NGC 6701 


SB(s)a 


57.2 


1.5 X 1.3 


-20.8 


0.01323 


3969 


32 


24 


Draco 


NGC 7770 


SO 


58.7 


0.8 X 0.7 


-19.4 


0.01414 


4242 


27 


50 


Pegasus 


NGC 7771 


SB(s)a 


60.8 


2.5 X 1.0 


-20.8 


0.01445 


4335 


66 


68 




Stephan's Quintet 
NGC 7317 


E4 


93.3 


1.1 X 1.1 


-20.3 


0.02201 


6603 


12 


150 


Pegasus 
jj 


NGC7318A ... 


E2 pec 


93.7 


0.9 X 0.9 


-20.5 


0.02211 


6633 








NGC7318B ... 


SB(s)bc pec 


82.0 


1.9 X 1.2 


-20.6 


0.01926 


5778 








NGC 7319 


SB(s)bc pec 


95.4 


1.7 X 1.3 


-20.8 


0.02251 


6753 


41 


148 




NGC 7320 


SA(s)d 


13.7 


2.2 X 1.1 


-17.5 


0.00262 


786 


59 


132 


5J 



six observing nights and 34 different pointings. The cen- 
tral position was observed in dithering mode to gain spatial 
resolution, while the remaining 33 positions were observed 
without dithering due to the large size of the mosaic. Seven 
positions were observed on the 28th October 2006, 19 posi- 
tions were observed between the 10th and 12th of Decem- 
ber 2007, 1 position on August 9th 2008 and the remaining 
pointings on October 28th 2008. Figure 1 shows the mosaic 
pattern covering NGC 628 consisting in a central position 
and consecutive hexagonal concentric rings. The area cov- 
ered by all the observed positions accounts approximately 
for 34 arcmin^, making this galaxy the largest area ever 
covered by a IFU mosaicking. The spectroscopic mosaic con- 
tains 11094 individual spectra, considering overlapping and 
repeated exposures (see Paper II). 



4.2 NGC 1058 

NGC 1058 is a well studied Sc spiral with a projected size 
of 3.0 X 2.8 arcmin at a distance of 10.6 Mpc, in the con- 
stellation of Perseus. The observations for this galaxy were 
performed on three consecutive nights from the 7th to the 
9th December 2007. The mosaic consists of the central po- 



sition and one concentric ring, covering most of the galaxy 
surface within one optical radius (defined by the B-band 
25th magnitude isophote). Ferguson et al. (1998) found the 
existence of H II regions out to and beyond two optical radii 
in this galaxy. We tried to observe these intrinsically inter- 
esting objects by performing a couple of offsets of 2 and 
2.5 arcmin north-east from the central position. These 2 ad- 
ditional position were merged to the original 7 tiles for a 
mosaic, covering an area of approximately 8.5 arcmin^ (see 
Figure 1). All positions (with the exception of one blind off- 
set) were observed in dithering mode, accounting for a spec- 
troscopic mosaic containing 7944 individual spectra. At the 
time of the observations, we were able to observe the recently 
discovered supernova 2007gr, a SN type Ic located at 24". 8 
west and 15". 8 north of the nucleus of NGC 1058 between 
two foreground stars (see subsection 7.4 and Figure 11). 



4.3 NGC 1637 

NGC 1637 is a SAB distorted galaxy in Eridanus with a 
projected size of 4.0 x 3.2 arcmin, at a distance of 12 Mpc. 
This galaxy presents a clear asymmetry with a third well- 
defined arm seen in optical images, an unusually extended 
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Figure 2. Examples of the spectroscopic mosaics of NGC 628, NGC 1058, NGC 1637 and NGC 3184. Eacli panel shows an intensity level 
narrow-band map centered at Ha (6563 A) in units of 10~^^ ergs"^ cm~^ arcsec"^. Note the effect of surface area coverage for dithered 
(e.g. NGC 1058, NGC 1637) vs. non-dithered (e.g. NGC 628, NGC 3184) observations, where the gaps between the fibres are clearly seen. 



HI envelope [Dhi /-D25 = 3.0), and an optical centre that 
differs from the kinematic centre by 9 arcsec (Roberts et al. 
2001). NGC 1637 was observed during December 8th to 10th 
2007. The mosaic was built with a central position and one 
concentric ring of 6 pointings (see Figure 1), covering most 
of the galaxy surface within one optical radius. The mosaic 
covers approximately 7 arcmin^ . This galaxy has a full spec- 
troscopic mosaic containing 6951 individual spectra. 

4.4 NGC 2976 

NGC 2976 is a SAc peculiar spiral galaxy with strong emis- 
sion line spectra with a projected size of 5.9 x 2.7 arcmin 
in Ursa Major, at a distance of 3.6 Mpc, being the clos- 
est object of the sample. The observations for NGC 2976 
were carried out on October 30th 2008. Given the distorted 
morphology of the galaxy a more convenient mosaic pattern 



was designed. Two pointings were observed for this object, 
corresponding to the central region of NGC 2976. The obser- 
vations were performed in non-dithering mode. The spectro- 
scopic data for this galaxy consist of 662 individual spectra. 

4.5 NGC 3184 

NGC 3184, a SAB face-on galaxy located in Ursa Major, 
has the 2nd largest angular size in the sample. It covers an 
area of 7.4 x 6.9 arcmin at a distance of 11 Mpc. NGC 3184 
has been classified as one of the metal-richest galaxies ever 
observed (McCall et al. 1985; van Zee et al. 1998), and 
which has also harboured recently a supernova explosion 
(SN 1999gi) (Nakano & Kushida 1999). Three concentric 
rings are necessary to cover the entire optical disk. Observa- 
tions for this galaxy were performed on December 10th 2007, 
following the standard mosaicking pattern with a central po- 
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Figure 3. Spectroscopic mosaics of NGC 3310, NGC 6701, NGC 7770, and the Steplian's Quintet. Note the eflfect of dithered (top panels) 
vs. non-dithcred (bottom panels) observations. Units as in Figure 2. 



sition and one complete ring of 7 IFU pointings. Then, on 
April 27th and 28th 2009, 9 additional pointings were ob- 
served covering partially a second concentric ring as shown 
in Figure 1. The area covered by all the observed positions 
is ~ 16 arcmin'^. The spectroscopic data for this galaxy con- 
sists of 5296 individual spectra. 



4.6 NGC 3310 

NGC 3310 is a very distorted spiral galaxy with strong star 
formation in the constellation of Ursa Major, at a distance 
of 17.5 Mpc. It covers an area of 3.1 x 2.4 arcmin in the op- 
tical -B-band, with a very bright central nucleus, surrounded 
by a ring of luminous HII regions. Different studies of this 
galaxy suggest a recent merging episode which triggered the 
burst of star formation (Kregel & Sancisi 2001; Wehner et al. 
2006, and references atherein). Given its morphology, a tai- 
lored mosaic pattern was constructed for this galaxy (see 



Figure 1). Three pointings cover the surface of NGC 3310 
with a central position centered in the galaxy's nucleus and 
two offsets of (-35, 35) and (35, -35) arcsec in (RA, Dec) in 
north-west and south-east directions respectively. The ob- 
servations were carried out on December 8th 2007, and were 
performed in dithering mode. This galaxy has a full spec- 
troscopic mosaic, which covers an area of approximately 2.8 
arcmin^. The spectroscopic data for this galaxy totals 2979 
individual spectra. 



4.7 NGC 4625 

NGC 4625 is a low-luminosity SAB, one-armed Magellanic 
spiral galaxy thought to be interacting with the also single- 
armed spiral NGC 4618 in Canes Venatici, at a distance of 
9 Mpc. The optical size of this galaxy covers an area of 
approximately 2.2 x 1.9 arcmin, however Gil de Paz et al. 
(2005) discovered an extended UV disk reaching to 4 times 
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Table 2. Summary of observations of the PINGS sample. Col. (1): 
Galaxy name. Col. (2): Number of individual IFU positions ob- 
served for each galaxy mosaic. Col. (3): Status of the mosaicking; 
when a percentage is shown, the number in parentheses repre- 
sents the total number of pointings necessary to cover the optical 
surface of the galaxy. Col. (4): Total number of spectra for each 
mosaic. Col. (5): Comments: a) Largest area ever covered by a IFU 
mosaicking, with a total of ~35 arcmin^; b) An offset position of 
NGC 1058 was not observed in dithering mode (see Figure 2); c) 
Galaxy with the second largest area of the sample; d) NGC 7770 
within the field; e) NGC 7318A and NGC 7318B are included in 
one field; *A11 pointings observed in dithering mode. 

Object Positions Mosaic Spectra Notes 



(1) (2) (3) (4) (5) 



NGC 628 


34 


92% (37) 


13571 


a 


NGC 1058 


9 


complete 


7944 


*,b 


NGC 1637 


7 


complete 


6951 


* 


NGC 2976 


2 


22% (9) 


662 




NGC 3184 


16 


84% (19) 


5296 


c 


NGC 3310 


3 


complete 


2979 




NGC 4625 


1 


14% (7) 


993 




NGC 5474 


6 


86% (7) 


5958 


* 


NGC 6643 


3 


complete 


2979 




NGC 6701 


1 


complete 


993 




NGC 7771 


3 


complete 


993 


d 


Stephan's Q. 


4 


complete 


1324 


e 



its optical radius showing evidence of recent star formation. 
The observation of this galaxy was performed on December 
9th 2007 with one single pointing in dithering mode covering 
the optical radius of NGC 4625. The spectroscopic data for 
this object consists of 993 individual spectra. 

4.8 NGC 5474 

NGC 5474 is a strongly lopsided spiral galaxy covering an 
area of 4.8 x 4.3 arcmin in Ursa Major, at a distance of 7 
Mpc. We observed this galaxy with a standard mosaic con- 
figuration consisting of one central position and one con- 
centric ring. All pointings were observed in dithering mode. 
Observations were carried out in two different periods; 4 
positions were observed during August 9th and 10th 2008, 
while 2 additional pointings were observed on the 27th April 
2009. Given the distorted morphology of this galaxy, the 
central position of the mosaic does not coincide with the 
bright bulge; a 30 arcsec offset in declination was performed 
towards the south, so that the area covered by the IFU mo- 
saicking includes most of the optical disk of the galaxy in 
a symmetric way. The area covered by all the observed po- 
sitions is approximately 6 arcmin^. The spectroscopic data 
for this galaxy totals 5958 individual spectra. 

4.9 NGC 6643 

NGC 6643 is a SAc galaxy in Draco, with a projected size 
of 3.8 X 1.9 arcmin in the B-band, at a distance of 20 Mpc. 
A tailored mosaic pattern for this galaxy was constructed in 
order to cover most of its optical area. Three pointings cover 
the surface of NGC 6643 with a central position centered on 
the bulge and two offsets of (37, 34) and (-35, -34) arcsec 



in (RA,Dec) in north-east and south-west directions respec- 
tively (see Figure 1). Observations were performed on June 
2nd 2008 for the first 2 positions and on August 10th 2008 
for the 3rd position, all of them in dithering mode. At the 
time of the first observing run, we were able to observe the 
supernova 2008bo, a SN type lb located at 31" north and 15" 
west of the nucleus of NGC 6643. This galaxy has a complete 
spectroscopic mosaic covering an area of approximately 2.8 
arcmin^. The data consists of 2979 individual spectra. How- 
ever, due to a technical problem with the instrument set-up, 
positions 1 and 2 do not cover the usual wavelength range, 
but are shifted towards the red by approximately 100 A. 

4.10 NGC 6701 

NGC 6701 is a small barred spiral in Draco with an angu- 
lar size of 1.5 X 1.3 arcmin, at a distance of 57 Mpc. This 
galaxy was considered to be an isolated galaxy, but studies 
of NGC 6701 have discovered morphological and kinemat- 
ical features that are consequence of an interaction, most 
probably with a companion at 73 kpc in projected distance 
(Marquez et al. 1996). The observation of NGC 6701 was 
carried out on August 9th 2008 with one single pointing in 
dithering mode covering the optical radius of the galaxy (see 
Figure 1). The spectroscopic data of this galaxy contains 993 
individual spectra. 

4.11 NGC 7770 and NGC 7771 

The main target for this mosaic was the galaxy NGC 7771, 
a barred spiral in Pegasus with an optical B-band size of 
2.5 X 1.0 arcmin at a distance of 59 Mpc. This galaxy is 
part of an interactive system containing mainly the face-on 
spiral NGC 7769 and the faint lenticular NGC 7770 (Nord- 
gren et al. 1997). The central part of NGC 7771 contains 
a massive circumnuclear starburst which was probably trig- 
gered by the interaction with the other members of the group 
(Smith et al. 1999, and references therein). Due to the pro- 
jected size of this galaxy, the mosaic pattern was constructed 
with three IFU positions. The central position of the mo- 
saic has an offset of (-15, -15) arcsec in (RA, Dec) from the 
geometrical centre of the galaxy (see Figure 1). Two ad- 
ditional positions were observed with offsets of (37, 33) and 
(-37, -33) arcsec. A second member of the interacting group, 
NGC 7770, a smaU SO galaxy with an optical size 0.8 x 0.7 
arcmin was observed within the field of the mosaic pattern. 
Observations of all 3 positions were performed on October 
30th 2008. The spectroscopic data for this galaxy contains 
993 individual spectra. 

4.12 Stephan's Quintet 

The Stephan's Quintet is one of the most famous and well- 
studied group of galaxies, consisting of NGC 7317, 7318A, 
7318B, 7319 and 7320 in Pegasus. The distance to this com- 
pact group of galaxies has been in debate due to the pres- 
ence of the brightest member, NGC 7320, which exhibits a 
smaller redshift than the others, suggesting that is a fore- 
ground object lying along the line of sight of the other four 
interacting galaxies. Although some controversy prevailed 
(Balkowski et al. 1974; Kent 1981), recent observations by 
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HST show that individual stars, clusters, and nebulae are 
clearly seen in NGC 7320 and not in any of the other galax- 
ies (Gallagher et al. 2001; Appleton et al. 2006, and refer- 
ences therein). Four individual pointings in non-dithering 
mode were observed on August 10th 2008, three of which 
were centered at the bright bulges of NGC 7317, 7319 and 
7320, while the last pointing was centered in configuration 
to cover NGC7318A and NGC7318B (see Figure 1). The 
spectroscopic data for all pointings of the Stephan's Quin- 
tet contains 1324 individual spectra. 



5 DATA REDUCTION 

The reduction of IFS observations possesses an intrinsic 
complexity given the nature of the data and the vast amount 
of information recorded in a single observation. This com- 
plexity is increased if one considers creating an IFU spectro- 
scopic mosaic of a given object for which the observations 
were performed not only on different nights, but even in 
different years, with dissimilar atmospheric conditions, and 
slightly differing instrument configurations. 

In this section we give an overview of the IFS data re- 
duction process for all the observations of the PINGS sam- 
ple. In general, the reduction process for the all pointings 
follows the standard steps for fibre-based integral field spec- 
troscopy. However, the construction of the mosaics out of 
the individual pointings requires further and more compli- 
cated reduction steps than for a single, standard IFU obser- 
vation. These extra steps arise due to the special mosaick- 
ing pattern for some of the objects, the differences in the 
atmospheric transparence and extinction, slight geometrical 
misalignments, sky-level variations, differential atmospheric 
refraction, etc. A complete explanation of the complex data 
processing for the creation of the PINGS mosaics sample 
is beyond the scope of this paper, but the reader will find 
a detailed description of the IFS data reducing in Sanchez 
2006 (hereafter San06) and additional information on the 
mosaicking technique for the PINGS sample in the descrip- 
tion of the PPAK-IFS survey of NGC 628 (Paper II). 

Following San06, all the data reduction steps can be 
summarised as follows: a) Pre-reduction. b) Identification of 
the location of the spectra on the detector, c) Extraction 
of each individual spectrum, d) Distortion correction of the 
extracted spectra, e) Application of wavelength solution, f) 
Fibre-to-fibre transmission correction, g) Flux calibration, 
h) Allocation of the spectra to the sky position, i) Cube 
and/or dithered reconstruction (if any). 

The raw data extracted from a PINGS observation con- 
sists of a collection of spectra, stored as 2D frames, aligned 
along the dispersion axis. The pre-reduction of the IFS data 
consists of all the corrections applied to the CCD that are 
common to the reduction of any CCD-based data, i.e. bias 
subtraction, flat fielding (in the case of PINGS, using twi- 
light sky exposures), combinations of different exposures 
of the same pointing and cosmic ray rejection. The pre- 
reduction processing was performed using standard IRAF^ 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



packages for CCD pre-reduction steps while the main re- 
duction was performed using the R3D software for fibre-fed 
and integral-field spectroscopy data (San06) in combination 
with the E3D visualization software (Sanchez 2004). 

On a raw data frame, each spectrum is spread over a 
certain number of pixels along the cross-dispersion axis. The 
spectra are generally not perfectly aligned along the disper- 
sion axis due to the configuration of the instrument, the 
optical distortions, the instrument focus and the mechanical 
flexures. Therefore, in order to flnd the location of each spec- 
trum at each wavelength along the CCD and to extract its 
corresponding flux, we made use of continuum illuminated 
exposures taken at each pointing corresponding to a different 
orientation of the telescope. Each spectrum was extracted 
from the science frames by co-adding the flux within an aper- 
ture of 5 pixels assuming a cut across the cross-dispersion 
axis found by iterative Gaussian fits (see subsection 6.5 for 
a detailed description of this reduction step). Since the mis- 
alignments of the fibres with the pseudo-slit also affect the 
wavelength solution, we require lamp exposures obtained at 
each observed position to find a wavelength solution for each 
individual spectrum. Wavelength calibration was performed 
using HeHgCd-|-ThAr arc lamps obtained through the in- 
strument calibration fibres. Differences in the fibre-to-fibre 
transmission throughput were corrected by creating a mas- 
ter fibrefiat from twilight skyfiat exposures taken in every 
run. 

The reduced IFS data can be stored using different data 
formats, all of which should allow to store the spectral in- 
formation in association with the 2D position on the sky. 
Two data formats are widely used in the IFS community: 
datacubes (3-dimensional images) and Row-Staked-Spectra 
(RSS) files. Datacubes are only valid to store reduced data 
from instruments that sample the sky-plane in a regular-grid 
or for interpolated data. In this case the data are stored in a 
3-dimensional FITS image, with two spatial dimensions and 
one corresponding to the dispersion axis. RSS format is a 
2D FITS image where the X and Y axes contain the spec- 
tral and spatial information respectively, regardless of their 
position in the sky. This format requires an additional file 
(either a FITS or ASCII table), where the position of the 
different spatial elements on the sky is stored. RSS is widely 
used by IFUs with a discontinuous sampling of the sky, as 
it is the case for PPAK. We chose to store the PINGS data 
in the RSS format, with corresponding position tables. 

Once the spectra are extracted, corrected for distor- 
tions, wavelength calibrated, and corrected for differences in 
transmission fibre-to-fibre, they must be sky-subtracted and 
flux-calibrated. One of the most difficult steps in the data 
reduction is the correct subtraction of the night sky emis- 
sion spectrum. In long-slit spectroscopy the sky is sampled 
in different regions of the slit and a median sky spectrum 
is obtained by spectral averaging or interpolation. This is 
possible due to the size of the long-slit compared with the 
size of the astronomical objects of interest. However, in IFS 
the techniques vary depending on the geometry of the ob- 
served object and on the variation of the sky level for a 
given observation. By construction, many of the positions 
of the PINGS mosaics (specially those in the galaxy centre) 
would fill the entire FOV of the IFU, and none of the spax- 



© 2010 RAS, MNRAS 000, 1-26 



PINGS: the PPAK IPS Nearby Galaxies Survey 11 



els^ would be completely free of galaxy "contamination" . In 
this case, we obtained supplementary sky exposures (imme- 
diately after the science frames) applying large offsets from 
the observed positions and between different exposures, we 
then used these "sky-frames" to perform a direct sky sub- 
traction of the reduced spectra. On the other hand, if the 
FOV is not entirely filled by the object, it is possible to select 
those spaxels (i.e. the sky-fibres in the case of PPAK) with 
spectra free of contamination from objects, average them 
and subtract the resulting sky-spectrum from the science 
spectra. We used this technique for observations in the last 
ring of a mosaic or at the edges of the optical surface of the 
galaxies, where the sky-fibres bundles did actually sample 
the sky emission. 

Once the sky emission has been subtracted, we need 
to flux calibrate the observed frames. Absolute spectropho- 
tometry with fibre-fed spectrographs is rather complex; as 
in slit-spectroscopy, where slit losses impose severe limita- 
tions, IFU spectrographs can suffer important light losses 
due to the geometry of the fibre-arrays. The flux calibra- 
tion requires the observation of spectrophotometric standard 
stars during the night. Given that the PPAK IFU bundle 
does not cover the entire FOV due to gaps among the fi- 
bres, the observation of calibration standard stars is prone 
to flux losses, especially when the standards are not com- 
pletely well centered in a single IFU spaxel. However we 
developed a method which takes into account the flux losses 
due to the gaps in the fibre-bundle, the pointing misalign- 
ments and PSF of the observed standard stars, as well as 
corrections for minor cross-talk effects, airmass, local opti- 
cal extinction and additional information provided by broad 
and narrow-band imaging photometry in order to obtain the 
most accurate possible spectrophotometric fiux-calibration 
within the limits imposed by the instrumentation. 

A total of six standard stars from the Oke spectropho- 
tometric candles (Oke 1990) were observed for the purpose 
of flux calibration during the observing runs. These frames 
were reduced following the basic procedure described above. 
To counteract the loss of flux, the observed spectrum of a 
standard star was obtained by adding up the spectra from 
consecutive concentric spaxel rings centered on the flbre 
where most of the standard's flux is found, until a conver- 
gence limit was found (see Figure 4). We then determined a 
night sensitivity curve as a function of wavelength by com- 
paring the observed flux with the calibrated spectrophoto- 
metric standard spectrum considering the flUing factor of the 
flbre-bundle. We applied this function to the science frames 
considering corrections for the airmass and the optical ex- 
tinction due to the atmosphere as a function of wavelength 
to get relative flux-calibrated spectra for each position. 

Figure 6 shows the shape and relative magnitude of nine 
sensitivity curves obtained during the three years of obser- 
vations. The difference in the vertical scale reflects the varia- 
tion of the spectrophotometric transmission during different 
nights and observing runs (the flux calibration obtained by 
applying these sensitivity curves is just a relative one, an 
absolute calibration is obtained by re-scaling by a factor de- 
rived after the comparison with the broad-band imaging by 
applying the method described below). However, the differ- 

^ Definition of IFS discrete spatial elements 
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Figure 4. Effect of standard star flux loss due to bad seeing and 
misalignments. On the top panel, the standard star BD-f 25d4655 
was observed during good atmospheric conditions and was well 
centred on the central fibre of the PPAK bundle. The figure con- 
tains 37 spectra corresponding to the 36 concentric fibres plus the 
central one. Most of the flux falls in the central flbres (red colour 
spectra in the online version), while very small residuals are seen 
in the rest of them (green- blue colour spectra). The bottom panel 
shows the spectra of the same star observed during turbulent at- 
mospheric conditions and when the star was not well centred in 
the flbre bundle. In this case the flux is spread over a large num- 
ber of concentric flbres and the individual contribution of each of 
them is important for the total observed flux. 



ences on the shape from one sensitivity curve to another as 
a function of wavelength, reflect the intrinsic dispersion of 
the flux calibration. The bottom panel of Figure 6 shows 
the variation of the sensitivity curves as a function of wave- 
length after a grey shift with respect to an arbitrary spline 
fltting normalised at the wavelength of H/3 (4861 A) . This 
normalisation wavelength was chosen as most spectroscopic 
studies normalise the observed emission line intensities to 
the flux in H/3. The maximum variation from the blue end 
of the spectrum compared to the red one is of the order of 
~ 0.15 mag, corresponding to a maximum calibration error 
of 15% due solely to the intrinsic dispersion in the rela- 
tive fiux calibration. However, the actual RMS is less than 
0.1 mag, corresponding to a typical error in the relative fiux 
calibration of less than 10%. 

For those galaxies with suitable multi-band photomet- 
ric data available in the literature (e.g. NGC 628, NGC 3184, 
NGC4625, NGC 5474), we used the first-order calibration of 
either the central pointing of a mosaic or the position with 
the highest S/N and best sky-subtraction to perform an ad- 
ditional correction to the absolute flux calibration by com- 
paring the convolved spectra observed in this held (taking 
into account fibre apertures and filters' response functions) 
to the corresponding flux measured hy B,V , R broad-band 
and Ho? narrow-band imaging photometry for the same po- 
sition. This procedure ensures a very precise flux calibra- 
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Figure 5. Observed standard stars on the 9th December 2007. 
The top panel figure shows a comparison of the night flux cal- 
ibration applied to the observed stars (in colour in the online 
version) and the corresponding calibrated fluxes. The panel be- 
low show the residuals in magnitudes as a function of wavelength 
for each star. The points with relatively large deviations are due 
to strong sky emission lines and cosmetic defects of the CCD. 



Figure 6. Variation of the spectrophotometric transmission as 
a function of wavelength. The top panel shows nine sensitivity 
transmission curves using different calibration standard stars. The 
thick dark line is a 3rd degree spline fitting to the data. The panel 
below shows the dispersion in magnitudes after a grey re-scaling of 
all the sensitivity curves with respect to the spline fit normalised 
at the H/3 wavelength (4861 A). 



tion and sky extinction correction for this master pointing. 
To our knowledge, no other IFS observations have ever at- 
tempted to get such (instrument-limited) spectrophotome- 
try accuracy. All galaxies belonging to the SINGS sample 
were corrected by this method using their ancillary data as 
described in more detail in Paper II. 

After reducing each individual pointing with a first- 
order flux calibration and with the help of the absolute flux- 
calibrated master pointing, we built a single RSS flle for the 
whole mosaic following an iterative procedure. The process 
starts in the master pointing chosen for a speciflc mosaic, 
i.e. the pointing that has the best possible flux calibration 
and sky extinction correction, with the best signal-to-noise 
and the most optimal observing conditions regardless of the 
geometric position of the pointing in the mosaic. Taking this 
master pointing as a reference, the mosaic is constructed by 
adding consecutive pointings following the particular mo- 
saic geometry. During this process, the new added pointing 
is re-scaled by using the average ratio of the brightest emis- 
sion lines found in the overlapping spectra (which is then 
replaced by the average between the previous pointing and 
the new re-scaled spectra). In most cases the scale factor is 
found to be between 0.7 and 1.3 with respect to the master 
pointing. 

However, this ratio is wavelength dependent (specially 
in the cases of variable photometric conditions between the 
observations). Therefore as a second-order correction, we fit- 
ted the variations found between the previous pointing and 
the new re-scaled overlapping spectra to a low order polyno- 
mial function and divided all the spectra in the new pointing 
by the resulting wavelength dependent scale. This correction 
has little effect on the data when the observations were per- 



formed during the same or consecutive nights, as it is the 
case for the small mosaics. However, we accounted variations 
after all the possible corrections of the order of 10-15% in 
the extreme cases when the observations were carried out 
at different epochs (e.g. NGC628, NGC3184). This level of 
error is what we expect from observations performed during 
different nights and observing runs, reflecting the variation 
of the spectrophotometric transmission from night to night 
(see Figure 6). The procedure was repeated for each mosaic 
until the last pointing is included (except for the Stephan's 
Quintet, where not actual overlapping occurs), ending with 
a final set of individual spectra and their corresponding po- 
sition tables. This process ensures a homogenous fiux cal- 
ibration and sky extinction correction for the entire data 
set. 



6 ERRORS AND UNCERTAINTIES IN THE 
DATA SAMPLE 

During the process of data reduction and basic analysis, we 
have identified several possible sources of errors and uncer- 
tainties in the PINGS data set. Each of them contribute 
in a different way and magnitude to the overall error bud- 
get associated with the observations. These are in order of 
importance: 1) sky subtraction, 2) flux calibration, 3) differ- 
ential atmospheric refraction (DAR), 4) cross-talk, and 5) 
second order spectra. In this section we describe the nature 
of each of these sources of errors, the tests performed in or- 
der to understand their effects on the accuracy of the data, 
and the techniques applied to minimise them. 
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6.1 Sky subtraction 

As mentioned before, sky subtraction is one of the most 
difficult steps in tlie IFS reduction process and it is partic- 
ularly complex for the nature of the observing technique of 
PINGS. As described in Sanchez (2006), a deficient sky sub- 
traction in this sort of data has several consequences: "the 
contamination of the sky emission lines along the spectra 
which prevents the detection and/or correct measurement 
of relatively weak nebular emission lines (e.g. the weak tem- 
perature sensitive [O III] A4363, which is located in the same 
spectral region as the strong Hg I A4358 sky line), and also 
affects the shape and intensity of the continuum, which is 
important for the analysis of the stellar populations and the 
determination of reddening". In fact, we made use of the 
mosaicking method in order to find the best possible sky 
subtraction. Due to the shape of the PPAK bundle and by 
construction of the mosaics in the standard mosaic configu- 
ration, II spectra of a given pointing (corresponding to one 
edge of the hexagon) overlap with the same number of spec- 
tra from the previous pointing (see NGC628 or NGC3184 
in Figure 2). This allows the comparison of the same ob- 
served regions at different times and with different atmo- 
spheric conditions. 

For a non-standard configuration the number of over- 
lapping fibres is larger (e.g. NGC3310). These overlapping 
spectra can be compared and used to correct for the sky 
emission of the adjacent frame. However, prior to perform- 
ing the sky subtraction it is required to visually check that 
no residual of the galaxy is kept in the derived spectrum. 
This can be the case if the transmission changed substan- 
tially during the observation of the adjacent frames. These 
techniques proved to result in good sky subtraction in most 
cases. On the other hand, when we were forced to obtain 
supplementary large-angle offset sky-exposures for the in- 
ner pointings in the mosaics that were completely filled by 
the target, we found that when the sky exposure is taken 
within a few minutes of the science exposure it produces a 
good subtraction. For those cases in which the atmospheric 
conditions changed drastically and/or the sky subtraction 
appeared to be poor, we combined different sky frames with 
different weights to derive a better result. 

One way of assessing the goodness of the sky subtrac- 
tion is to check for sky residuals in the subtracted spectra. 
The galaxy mosaic more prone to be affected by residuals 
in the sky subtraction is NGC 628, which as explained in 
section 4, was observed during six nights along four observ- 
ing runs. Therefore we would expect that the spectroscopic 
mosaic of this galaxy would show the most extreme effects 
due to the sky subtraction to be found in the PINGS sam- 
ple, given all the variations in transparency and photometric 
conditions of the night-sky along the three years of observa- 
tions. 

In order to obtain a quantitative assessment of the qual- 
ity of the sky subtraction, we performed two different data 
reductions of the spectroscopic mosaic of NGC 628. In the 
first reduction, the sky subtraction was performed directly 
with the average spectrum of the sky fibres at each posi- 
tion, without considering the overlapping spectra between 
pointings, and not accounting for the object contamination 
in the sky fibres. Therefore, in this first reduction we applied 
a "poor" sky subtraction. For the second reduction, we ap- 



plied an individual sky subtraction per mosaic position using 
the techniques described above, i.e. applying corrections us- 
ing the overlapping spectra, checking for galaxy residuals in 
the derived sky spectrum, using the sky exposures obtained 
by large offsets for the most internal regions of the galaxy, 
and combining different sky frames with different weights in 
those cases when there were important changes in the trans- 
mission between pointings observed during the same night. 
We refer to this reduction as the "refined" sky subtraction. 

Airglow is the most important component of the light of 
the night-sky spectrum at Calar Alto observatory, although 
a substantial fraction of the spectral features is due to air 
pollution (Sanchez et al. 2007). The strongest sky line in 
the Calar Alto night-sky spectrum is the [O I] A5577 line, 
followed by the [O I] A6300 line, both produced by airglow 
with a notorious stronger effect near twilight. A deficient 
sky subtraction can be recognised by residual features of 
the sky lines in the derived spectra, this effect is clearly 
seen in the [O I] A5577 sky line which is located in a spec- 
tral region without any important nebular emission line. In 
general terms, (without considering variations in the trans- 
parency of the sky) , a residual in emission of this line would 
imply a subtraction of the sky spectrum of slightly lower 
strength than required, while an absorption feature would 
imply an over-correction. 

In order to make a comparative analysis of the strength 
of the sky residuals in the two data reductions of NGC 628 
described above, we measured the equivalent width (EW) 
of the residual features centered at the [O I] A5577 line. For 
numerical reasons (regions of null continuum), the local con- 
tinuum in the neighbourhood of the [O I] A5577 sky line was 
re-scaled to a flux level of 10~^® ergs~^ cm~^ A"'^ in every 
single spectrum of both mosaics. EWs with negative sign 
correspond to residual emission features, while positive EWs 
correspond to absorption features. 

Figure 7 shows the value of the EW residuals for the 
[O I] A5577 line for both data reductions as a function of the 
pointing position in the spectroscopic mosaic. Each posi- 
tion bin contains 331 values corresponding to the number of 
spectra per pointing, a total of 11104 values are shown, cor- 
responding to the 34 positions observed for NGC 628. The 
green dots correspond to the poor sky subtraction reduction, 
while the blue dots correspond to the refined sky subtraction 
(see the online version of this plot). There is a considerable 
amount of scatter of the EW residual value for the poor sky 
subtraction compared to the refined one. In the first two 
pointings of the mosaic (which correspond to central posi- 
tions of the galaxy), there are strong residuals in emission 
for the poor reduction, while the residuals have been min- 
imised in the refined one. Figure 7 shows clear evidence of 
those pointings in which the sky transparency varied by a 
considerable amount (positions 4, 5, 11, 12, 23, 24, 28, 29, 
33, 34). In all the pointings, the scatter in the residuals is 
improved in the refined reduction with respect to the first 
one. This effect is more notorious between positions 13 to 
22. The poor sky subtraction yields very strong sky residu- 
als in emission for positions 33 and 34, while in the refined 
reduction these are minimised. 

At the chosen continuum level used for this exercise, 
a (absolute) value of 5 A in EW for the [O I] A5577 resid- 
ual line in emission corresponds approximately to a ffux in- 
tensity value of 5 X 10~^^ ergs~^ cm~^, while a value of 
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Figure 7. Equivalent widtli values for the emission and absorp- 
tion residual features of the [O l] A5577 sky line as a function of 
the observed position for the spectroscopic mosaic of NGC 628. 
Positive values correspond to absorption features, while negative 
values to emission residuals. The local continuum level was re- 
scaled to the same value in order to make this comparison. The 
green dots (in the online version) correspond to a poor sky sub- 
traction, while the blue dots represent the refined reduction as 
explained in the text. The two horizontal dotted lines mark the 
threshold EW values of residuals features corresponding to a good 
sky subtraction. 

10 A corresponds to ~ 10 x 10~^^ ergs~^ cm~^. The aver- 
age flux intensity of the [Ol] A5577 sky line in Calar Alto 
is of the order of 33 x 10~^® ergs~^ cm"^ (Sanchez et al. 
2007). However, from a sample of 500 sky spectra acquired 
during the three years of observation we measured the in- 
tensity of the [Ol] A5577 in the range between 30 and 60 
X 10~^^ ergs~^ cm~^, with a mean value of 44. Therefore, a 
value of 5 A in EW for a residual emission feature would cor- 
respond to ~ 8-10% of the total emission of the [O I] A5577 
line. Visual inspection of the spectra with emission residual 
of the order of 5 A in EW confirms that this value could 
be considered as the threshold for a good sky subtraction. 
Spectra with emission or absorption residuals with absolute 
EW values less than 5 A could be considered to have a good 
subtraction, for features above this value the effects of a 
deficient sky subtraction are evident. 

The two horizontal dotted lines in Figure 7 indicate the 
±5 EW threshold value for both emission and absorption 
features. These two lines encompass a region for which the 
spectra can be considered with a good sky subtraction. The 
poor sky subtraction (green) shows a lot of scatter and a 
small fraction of the spectra falls within these limits. On the 
other hand, for the refined reduction (blue) a total of 9629 
spectra fall within these limits, i.e. 87% of the total mosaic. 
The number of spectra with sky subtraction problems for 
which |EW| > 5 A is 1475, i.e 13% of the mosaic, these spec- 
tra are found in those pointings with extreme transparency 
variations as expected. 

Figure 8 shows the histograms of the EW values for both 
data reductions, the poor sky subtraction in green and the 
refined reduction in blue colour following Figure 7 (see the 
online version of this plot). The top panels show the distri- 
bution of residual emission values, while the bottom panels 
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Figure 8. Histograms of the EW residual features shown in 
Figure 7. The top panels shows the distribution of EW values for 
the emission residuals. The bottom panels shows the EW values 
for the absorption features. All EW values are shown in logarith- 
mic scale (assuming a positive EW value for the omission resid- 
uals). The top histograms in both cases correspond to the poor 
sky subtraction, while the bottom histograms correspond to the 
refined reduction (green and blue colour histograms in the online 
version). The ~5 A EW threshold value is shown as the vertical 
line in all the histograms, residual values to the right of this line 
can be considered a deficient sky subtraction as explained in the 
text. 



show the absorption residual features for the [O I] A5577. 
The 5 A EW threshold value is shown as the vertical line 
in the histograms, residual values to the right of this line can 
be considered a deficient sky subtraction. Visual inspection 
of the spectra shows that, at the continuum level used for 
this comparison, emission or absorption features with values 
of log(jEW|) ^ could be considered negligible and within 
the statistical noise of the spectra. The residual emission 
histograms show that the poor sky subtraction produces a 
large number of strong residuals with values of |EW| > 5 
A, even reaching ]EW] ~ 60 A. The majority of the residual 
values in refined sky subtraction are found at log(|EW|) ~ 0, 
corresponding to negligible residual values, however there is 
a small tail of strong emission residuals for which |EW| > 5 
A (~ 18% of the total emission residuals). The distribution 
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of EW values of the absorption features for the poor subtrac- 
tion is approximately centered at the threshold limit, while 
for the refined reduction, the values are nearly normally dis- 
tributed with a centre value of log(EW) ~ with a small 
tail of strong absorption values 7%) due most likely to 
an over subtraction of the sky spectrum. 

The refined sky subtraction was the final adopted one 
for the spectroscopic mosaic of NGC 628. All the sky sub- 
traction techniques implemented showed that the quality of 
the derived spectra was improved by a considerable amount 
compared to a standard sky subtraction. Most of the sky 
residuals are within the limits of a reasonably good sky 
subtraction. The spectra with strong features are found for 
those positions in which the photometric conditions changed 
drastically during the night or observing run. This residual 
analysis allows to identify those pointings with strong sky 
variations and thus, to flag the spectral data for future anal- 
ysis. The sky subtraction for the rest of the PINGS sample 
was performed similarly to the refined technique described 
above. Therefore we applied the best possible sky subtrac- 
tion to all the spectroscopic mosaics within the limitations 
imposed by the IFS data itself. 



6.2 Detection of the [O III] A4363 line 

The Hg A4358 sky line strongly affects any attempt to mea- 
sure precisely the emission of the faint temperature-sensitive 
[OIII] A4363 line in any object with a low redshift. The fact 
that the strength of this line decreases with increasing abun- 
dance (Bresolin 2006), in combination with typically faint 
If II regions and low spectroscopic resolution limits impor- 
tantly the detectability of this key diagnostic line. 

In order to assess the significance of the detection of 
the [OIII] A4363 given the contamination of the Hg A4358 
sky line in our data, we performed a simulation of the de- 
tectability of the [O III] A4363 line for a given range of red- 
shifts and line intensity strengths. We simulated a pure emis- 
sion line spectrum including the H7 A4340 and [O III] A4363 
lines at the same spectral resolution of the PINGS observa- 
tions; we assumed a normally distributed I(A4363)/H7 ra- 
tio with a mean value of 0.10 ±0.05, corresponding to typ- 
ical values found in previous spectroscopic studies where 
the [Olll] A4363 line was detected in HII regions within 
the metallicity range of the PINGS sample (e.g. McCall 
et al. 1985); we did not consider higher ratios (~ 0.25 ± 0.10) 
which are representative of extremely low metallicity objects 
(e.g. Pagel et al. 1992; Izotov & Thuan 2004). We added a 
random statistical noise of 0.05 RMS at the continuum level 
constructed from the observed spectroscopic data. A sam- 
ple of 540 sky spectra were selected among all the observ- 
ing runs during the three years of observations (considering 
very different photometric conditions). A fiux calibrated sky 
spectrum was created out of these selected spectra. This sky 
spectrum was added to the previous emission line plus the 
noise described above to create a simulated "observed" spec- 
trum. An average sky spectrum constructed from a random 
subset of 36 sky spectra (the number of PPAK sky-fibres) 
was then subtracted from the simulated "observed" spec- 
trum to obtain a "sky-free" spectrum. 

Emission line intensities were then measured simultane- 
ously for both lines in the simulated "sky-free" spectrum us- 
ing the techniques described in subsection 7.2. These line in- 
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Figure 9. Detectability of the [OIII] A4363 line as a function of 
redshift and line strength according to the simulation described 
in subsection 6.2. For a given redshift, the [O III] A4363 can be 
significantly detected for flux values above the thick line, which 
marks the region at which the difference between the observed and 
simulated line intensity is of the order of 15%. The contamination 
effect of the Hg A4358 disappears for redshift values larger than 
0.004. 



tensities were then compared with the flux of the pure emis- 
sion lines. For a given redshift, we varied the emission line 
strengths of the simulated spectrum from high to lower val- 
ues until the significance of the detection of the [O III] A4363 
fell drastically. We performed 500 realisations of the emis- 
sion line intensity measurements for a given redshift and for 
a given line strength. Figure 9 shows the results of the simu- 
lation, the thick line represents the region at which the differ- 
ence between the line intensity measured from the simulated 
"sky-free" spectrum and the fiux from the pure emission line 
is of the order of 15%. According to the simulation, observed 
fiux values of the [O III] A4363 above this line can be signif- 
icantly detected at a given redshift. For fiux values below 
this region the significance of the detection is negligible as 
it is mostly affected by the statistical noise of the data. The 
contamination effect of the Hg A4358 disappears for redshift 
values larger than ~ 0.004, where the detectability of the 
[O III] A4363 depends on the signal-to-noise of the spectrum 
at low line-intensity levels. Experience with the data has 
proven that the simulation described in this section places 
very good limits on the detectability and potential measure- 
ment of this faint line, although after detection, individual 
and visual inspection of the spectra has to follow in order 
to correctly assess the usefulness of this line. 



6.3 Flux calibration 

Several refinements in the observation technique and stan- 
dard fiux recovery were applied to the pipeline which im- 
proved substantially the accuracy of the sensitivity functions 
obtained after every standard candle observation. During 
most of the observing runs we observed different standard 
stars per night at different airmasses in order to asses the 



© 2010 RAS, MNRAS 000, 1-26 



16 Rosales- Ortega et al. 



variation in transmission and its effect on the relative flux 
calibration. 

We generated several sensitivity curves following the 
standard pipeline procedure in R3D, changing the key pa- 
rameters that could affect the accuracy of the derived sen- 
sitivity function (e.g. order and type of the fitting function, 
extinction, airmass, smoothing, etc.). Furthermore, we made 
a comparison of the response curves obtained through R3D 
and the ones obtained using standard long-slit flux calibra- 
tion routines in IRAF after performing all the appropriate 
corrections and transformations for the two different kinds 
of data. We even derived whole-run sensitivity curves af- 
ter the combination of several response curves for a given 
observing run and applied the derived calibration to the ob- 
served standard candles as a proof of self-consistency. In all 
cases, we found very consistent results in the final relative 
flux calibration. As described in section 5, even without a 
re-calibration using broad-band imaging, the spectral shape 
and features are reproduced within the expected errors for 
an IFS observation (~ 20% in the absolute sense) along the 
whole spectral range, with a small increase in the blue region 
(A < 3800 A) due in part to the degradation of the CCD 
image quality and instrumental low sensitivity towards the 
blue (~ 2-5%, telescope/atmosphere excluded) in this spec- 
tral region (Roth, priv. comm.). 

For those galaxies with available multi-band photomet- 
ric data, small differences in the transmission curves of the 
filters and astrometric errors of the built mosaics bring some 
uncertainties in the derived flux ratios that contribute to the 
overall standard deviation when applying the photometric 
re-calibration described before. The errors in the flrst case 
are difflcult to estimate, however the latter ones were esti- 
mated by simulating different mosaic patterns by applying 
normally distributed random offsets of the simulated fibre- 
apertures with mean values of 0.3, 0.5 and 1.3 arcsec on 
the broad-band images and then comparing the extracted 
spectrophotometry. From the results of these simulations 
and considering that the IFS mosaics were re-centered us- 
ing the information directly from the aperture photometry, 
we expect that the location of the fibres lies within 0.5 arc- 
sec, and therefore the error due to the uncertainty in the 
astrometry would be of the order of 10%. Based on these 
results, we estimate a spectrophotometry accuracy better 
than ~ 0.2 mag, down to a flux limit corresponding to a sur- 
face brightness of ~22 mag/arcsec^ (see Paper II) when we 
apply the re-calibration derived by the flux ratio analysis. 
In the following section we compare our data with previ- 
ously published spectrophotometrically calibrated data. The 
spectral shape, the spectral features and emission line in- 
tensities match remarkably well, even for those objects for 
which no re-calibration was performed (see Figure 11). As 
stated above, to our knowledge no other IFU observation 
has ever implemented such corrections in order to obtain an 
instrumental-limited spectrophotometry accuracy as in this 
work. 

6.4 Differential Atmospheric Refraction 

An important systematic effect in any spectroscopic obser- 
vation is due to the refraction induced by the atmosphere, 
which tends to alter the apparent position of the sources 
observed at different wavelengths. By definition, there is no 



refraction when the telescope is pointed at the zenith, but for 
larger zenith angles the effect becomes increasingly signifi- 
cant. For IFU observations, this has the consequence that, 
when comparing for example the intensities at two different 
wavelengths (e.g the emission line ratio of a source), one will 
actually compare different regions, given that different wave- 
lengths are shifted relative to each other on the surface of 
the IFU. In theory, one is capable of performing a correction 
of DAR for a given pointing without requiring knowledge of 
the original orientation of the instrument and without the 
need of a compensator, as explained by Arribas et al. (1999). 
The correction of DAR is important for the proper combi- 
nation of different IFS exposures of the same object taken 
at different altitudes and under different atmospheric condi- 
tions, and for the proper alignment of a mosaic and dithered 
exposures, as it is the case of most PINCS observations. An 
IFS observation can be understood as a set of narrow-band 
images with a band-width equal to the spectral resolution 
(San06). These images can be recentered using the theoret- 
ical offsets determined by the DAR formulae (Filippenko 
1982) by tracing the intensity peak of a reference object in 
the FOV along the spectral range, and recentering it. The 
application of this method is basically unfeasible in slit spec- 
troscopy, which represents one additional advantage of IFS. 
The correction can be applied by determining the centroid of 
a particular object or source in the image slice extracted at 
each wavelength from an interpolated data cube. Then, it is 
possible to shift the full data cube to a common reference by 
resampling and shifting each image slice at each wavelength 
(using an interpolation scheme), and storing the result in 
a new data cube. A pitfall of this methodology is that the 
DAR correction imposes always an interpolation in the spa- 
tial direction as described above, a 3-dimensional (3D) data 
cube has to be created for each observed position, reducing 
the versatility given by the much simpler and handy RSS 
files. 

We have to note here that the widely accepted formula- 
tion summarised by the work of Filippenko and the concept 
of parallactic angle are just a first order approximation to 
the problem. All this theoretical body is based on the as- 
sumption that all different atmospheric layers have an equal 
refraction index, are flat-parallel and perpendicular to the 
zenith. While this approximation is roughly valid, there are 
appreciable deviations due to the topography and landforms 
at the location of the telescope, since they alter considerably 
the structure of the low-altitude atmospheric layers. There- 
fore, the "a posteriori" correction of the DAR effect, only 
possible when using IFS, is the most accurate approach to 
the problem. 

In general, the effects due to DAR in IFS are only im- 
portant for IFUs with small spatial elements (^J 1.5 arcsec) 
while for large ones (as it is the case of PPAK), the effect is 
mostly negligible, as experience with the instruments shows, 
especially when the airmass of the observations is 1.1 or be- 
low (Sandin et al. 2008). According to the DAR formulas one 
can calculate the angular separation in arcsec due to this ef- 
fect for two different wavelengths under typical atmospheric 
conditions for a range of airmasses (Filippenko 1982, Ta- 
ble 1). If we consider for example the wavelengths of Hq at 
A6563 and [O II] A3727 emission lines, the angular separation 
due to DAR is smaller than the radius of the PPAK flbres 
(1.35 arcsec) for airmasses below 1.3. Nevertheless, for each 
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object in the sample we analysed those pointings which were 
observed at an airmass > 1.2 in order to test any effects due 
to DAR in our data. We transformed these individual point- 
ings into 3D data cubes (with a scale of 1" /pixel) and we se- 
lected suitable sources within the field (e.g. foreground stars, 
compact bright emission line regions) to perform a DAR 
correction creating continuum maps of these bright sources 
and looking for spatial deviations along the dispersion axis. 
No significant intensity gradients were found in any of the 
test fields. Additionally, we looked for regions in which we 
could observe emission in the blue (e.g. [O II] A3727) but 
no emission in the red (e.g. Ha A6563) and vice versa. In 
most of the cases we did not find strange [Oil] /Ha ratios, 
although for some pointings we did not have enough bright 
HII regions to perform this exercise. However, for a number 
of pointings we did find peculiar deviations from the flux 
measured in both part of the spectra. All these pointings 
were observed with an airmass > 1.4 and correspond to: 
NGC 1637 (all 7 dithered pointings), NGC6701 (all 3 point- 
ings), and NGC 5474 (pointings 3, 4 & 5). These particular 
pointings could be individually converted to 3D data cubes 
and be further analysed and corrected for any DAR effects. 
In the case of particular studies which require a degree of 
spatial precision, one has to bear in mind the significance 
of the spatial-spectral information derived from the individ- 
ual fibres of these pointings. Integrated spectra over a large 
aperture (~ 5 arcsec) could be considered reliable. 

Considering in general that the effects of DAR are very 
small given the relatively large size of the spatial elements of 
PPAK, that ~ 70% of the observations were performed with 
an airmass 1.35, and that we found no significant evi- 
dence of DAR effects in the analysis of our data (apart from 
the flagged pointings described above), we decided that no 
DAR correction should be applied when building the flnal 
spectroscopic mosaics of the PINGS sample, avoiding the 
transformation of the RSS files to 3D data cubes and the 
undesirable interpolation of data. It is important to note 
that problems caused by atmospheric dispersion cannot be 
completely avoided in the spectroscopic study of extended 
objects like any other IPS observation, only for those ob- 
servations performed in dithering mode a correction for the 
effect of DAR can be sought. 

6.5 Cross-talk 

One of the critical reduction steps is to extract the flux corre- 
sponding to different spectra at each pixel along the disper- 
sion axis (see section 5). However, this is not straightforward 
as there could be contamination by flux coming from the ad- 
jacent spectra, i.e. cross-talk. This contamination effect may 
produce a wrong interpretation of the data. For observations 
performed with PPAK, and due to the geometrical construc- 
tion of the instrument, the adjacent spectra on the plane of 
the CCD may not correspond to nearby locations in the sky 
plane (Kelz et al. 2006). Furthermore, the position of the cal- 
ibration fibres (which are located in between the science ones 
along the pseudo-slit) also contributes to overall contamina- 
tion. Therefore, the cross-talk effect would potentially mix 
up spectra from locations that may not be physically related 
(because of the spectra position on the CCD) or from spec- 
tra with completely different nature (i.e. calibration fibres). 
Given that the cross-talk is an incoherent contamination it 
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Figure 10. Example of a strong cross-talk effect in a fibre ad- 
jacent to a very bright source, in this case a standard star with 
a maximum of 40000 counts. The top spectrum shows the full 
cross-talk effect when no correction is applied. The middle spec- 
trum shows a first-order correction obtained by minimising the 
pixel extraction aperture (but losing flux in the process). The 
bottom spectrum shows the correction of the cross-talk by the 
improved Gaussian-suppression method (described in Paper II) 
with an extracted spectrum of the same order of intensity as the 
average background signal (dotted line) found in the nearby fi- 
bres for this particular observation. The apparent periodicity of 
the cross-talk with respect to the wavelength is just an effect of 
the spectra misalignments and distortions on the CCD. 

is preferable to keep it as low as possible, an average value 
of ~ 1% with a maximum of ~ 10% (San06, and references 
therein) seems technically as a good trade-off. 

Given the limited size of CCDs and the need to record as 
many spectra as possible within that reduced area, IFS ob- 
servations will always face a certain level of cross-talk. The 
FWHM of the projected profile along the cross-dispersion 
axis is normally defined by the design of the spectrograph 
and the size of the input fibres, placing a limit in the se- 
lected aperture. As explained by San06, the PPAK spec- 
tra profiles have a FWHM of the order of ~ 2.3 pixels, and 
A peaks of ~ 5 pixels in the 2x2 binning mode (as was the 
case for all PINGS observations). Selecting an aperture size 
of the order of A peaks seems to be an acceptable compro- 
mise between maximising the recovered flux and minimising 
the cross-talk. Several methods have been implemented to 
minimise the effect of cross-talk, in particular, SanQ6 de- 
veloped an efficient technique named Gaussian-suppression 
that reduces the effects of the cross-talk and maximises the 
recovered fiux to within 10% of the original values for any 
spaxel at any wavelength. However, for certain raw frames 
which were too crowded or when we targeted bright sources 
within the field (e.g. foreground stars, galaxy bulges), we 
still found some level of contamination that could not be 
considered negligible. 

Therefore, we improved the Gaussian-suppression tech- 
nique to a new method which increases the signal-to-noise 
ratio of the extracted spectra and reduces the effect of 
the cross-talk compared to previous extractions. The new 
technique is explained in detail in Paper II, and it is now 
fully implemented in the standard R3D package (San06). 
We performed several tests using simulated and real spectra 
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with a broad range of intensities in order to assess the 
level of contamination that the extracted spectra show due 
to the cross-talk. We tested the new method varying the 
relative intensity of the spectra in the central and adjacent 
fibres, the pixel extraction apertures and the average width 
of the Gaussian profiles. We found that in the extreme cases 
(very bright adjacent fibres) the cross-talk is suppressed 
by 95%, being almost negligible for the range of spectral 
intensities compared to the instrumental misalignments 
and distortions found during conventional observations of 
the PINGS sample (see Figure 10). 

In summary, the two most important sources of error in 
the data reduction arise from the sky subtraction and the 
flux calibration. We have made several tests and applied new 
techniques in each case in order to minimise the magnitude 
of the uncertainties. We found that the sky subtraction can- 
not be applied as part of any standard reduction pipeline 
and has to be considered on an individual basis, depend- 
ing on the observing mode, configuration and strategy. On 
the other hand, the maximum expected error in the abso- 
lute flux calibration for those objects in which a broad-band 
imaging recalibration was applied is of the order of 20%, be- 
ing slightly larger in a narrow blue spectral region A < 3800 
A. For those objects without broad-band imaging, the abso- 
lute flux calibration error is of the order of 30%. However, 
the absolute error is better than 20% for objects observed 
during photometric conditions and/or with high S/N. The 
colors, spectral features and gradients are completely repro- 
duced when comparing our data with previously published 
long-slit observations (see subsection 7.1). 

DAR effects were found to be negligible due to the large 
size of the spatial elements of PPAK, we did not find strong 
evidence of DAR in the several test pointings analysed. How- 
ever we fiagged all individual IFU positions for which the 
observations were performed with an airmass > 1.3. With 
respect to cross-talk, the new extraction method proved to 
suppress this effect to a negligible level for the range of spec- 
tral intensities in the PINGS observations. Furthermore, no 
evidence of any second-order contamination has been found 
during the data reduction and analysis of the PINGS data 
set. As explained in section 4, the contamination of the 2nd 
order, up to 7200 A, is expected to be lower than 1/10000, 
being negligible for our science case. 



7 THE PINGS DATA 

The PINGS data set contains more than 50000 flux cali- 
brated spectra for a sample of 17 galaxies covering in total 
an observed spectroscopic area of nearly 80 arcmin^. Ta- 
ble 2 shows a summary of the observations, including the 
number of individual IFU pointings observed for the mosaic 
of each galaxy, the observational status of the mosaicking, 
the total number of spectra, and additional individual com- 
ments on each object. The spectroscopic data set samples 
the observed objects with flbre circular apertures of ^ 2.7 
arcsec in diameter, covering the optical wavelength range 
between ~ 3700 -7100 A, which includes the most promi- 
nent recombination and coUisionally excited emission lines 
from [Oil] AA 3727, 3729 to [SII] AA 6717, 6731. 

It is beyond the scope of this article to extract all the 



information potentially available in this huge spectroscopic 
data base. Independent studies for individual objects, re- 
gions and galaxy subsamples will be presented in a series of 
future papers. We consider here a few relevant scientiflc ex- 
amples that have been extracted from the data set to demon- 
strate the use of these data. In this section we present the 
integrated spectra and emission line intensities for selected 
H II regions of the IFS mosaics of NGC 1058 and NGC 3310, 
with a comparison to previously published data. Further- 
more, we obtained emission line maps for NGC 1058 and 
we present a qualitative description of the 2D distribution 
of physical properties inferred from the line intensity maps. 
Finally, we present the spectra of the SN 2007gr observed 
in NGC 1058 as an example of the utility of wide-fleld IFS 
surveys to find potential SNe progenitors and their environ- 
mental properties. 

7.1 Integrated spectra 

The study of the optical integrated light provides a pow- 
erful tool to investigate the physical properties of galaxies 
at different epochs in the history of the universe. Spectral 
investigations based on integrated optical spectrophotome- 
try have been used to explore the main drivers of galaxy 
evolution, e.g. star formation rate, star formation history, 
stellar mass, chemical abundance, dust content, etc. (e.g. 
Kauffmann et al. 2003; Tremonti et al. 2004). The analysis 
of the integrated spectra in nearby objects can be used to as- 
sess the limitations imposed by high-redshift surveys, such 
as their limiting magnitude and incomplete spatial cover- 
age (or aperture bias), factors that may be important given 
that many physical properties of galaxies vary depending on 
the geometry and position (e.g. stellar populations, metal- 
licity, extinction, etc.). A solution to the incomplete spa- 
tial coverage of spectroscopic observations on nearby ob- 
jects consists on using sequential methods that use time to 
scan a target while recording the spectral information. The 
standard method for obtaining integrated spectra in nearby 
objects was developed by Kennicutt (1992). The so-called 
drift-scanning technique consists of a scan perpendicular to 
the slit over the optical extent of the galaxy during a sin- 
gle exposure. Using this method Moustakas & Kennicutt 
2006b (hereafter MK06) obtained spatially integrated opti- 
cal (3600 - 6900 A) spectrophotometry for 417 nearby galax- 
ies of a diverse range of galaxy types, which were later used 
to study several integrated galactic properties (SFR, metal- 
licity, etc.). 

One by-product of 2D spectroscopy and IFS data sets 
is the intrinsic capability of adding up all the spectra within 
an observed field or mosaic into a single spectrum, i.e. us- 
ing the IFU as a large-aperture spectrograph to obtain the 
integrated spectra of a given FOV. The PINGS sample is 
an ideal data set for this purpose given that the spectro- 
scopic mosaics cover, in most cases, the entire optical radius 

* The analysis of this section was partly based on a set of 
IDL routines created in order to handle and visualise the 
2D spectroscopic mosaics. Although they were first devel- 
oped for the specific PINGS data format, they can be imple- 
mented to any IFS data based on RSS files and correspond- 
ing position tables. The PINGS software is freely available at 
http : //www. ast . cam. ac .uk/research/pings/. 
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Figure 11. Top panels: Integrated spectra of NGC 1058 derived from the PINGS slit-extraction (thick solid line) and the whole mosaic 
(dotted line) spectra compared to the MK06 data (blue line in the online version), as explained in the text. Flux in units of lO"'^^ 
erg s~^ cm~^ A~^. Bottom panels: Comparison of the integrated spectra of NGC 3310 with the same color coding as the figure above. 
For both galaxies, the diagrams on the right represent narrow-band intensity maps centred at Ha showing the mosaic pattern and the 
rectangular slit-apertures of MK06. The arrow in the NGC 1058 mosaic diagram shows the position of SN 2007gr. The contours of the 
individual fibres have been drawn in the case of NGC 3310. 



of the galaxy. The integrated spectra derived from PINGS 
can be used to study the real average spectroscopic proper- 
ties of a given nearby, large angular size galaxy, as opposed 
to previous studies that attempted to describe their aver- 
age properties by the analysis of individual spectra taken 
from different regions, or by targeting objects with a lim- 
ited extraction aperture which recovers only a fraction of 
the total optical light. In this section we present examples 
of high signal-to-noise integrated spectra for NGC 1058 and 
NGC 3310 obtained by co-adding the spectra from their cor- 
responding mosaics using different simulated apertures ap- 
plied to the IFS mosaic. We compare these data with previ- 
ously published integrated spectra from MK06. In Paper II, 
we present the integrated spectra of NGC 628 with a more 
elaborated analysis of the integrated stellar populations and 
nebular abundances. 

Figure 11 shows in the left upper panel the compari- 



son of the integrated spectrum of NGC 1058 extracted from 
the PINGS mosaic to the drift-scan spectrum obtained by 
MK06 of this galaxy, at the same spectral resolution. The 
right upper panel shows the PINGS spectroscopic mosaic 
diagram of NGC 1058 with intensity levels corresponding to 
a "narrow-band" 100 A width image extracted from the IFS 
data centred at Ha. The rectangular shape shows the slit- 
aperture used in the drift-scan technique by MK06 to obtain 
the integrated spectrum of this galaxy^. Three foreground 
stars and the supernova SN 2007gr (see subsection 7.4) were 
removed from the co-added spectra (the arrow shows the po- 
sition of SN 2007gr between two foreground stars). The black 



^ The simulated extraction slit was registered over the surface of 
the galaxy by visual comparison to corresponding diagram shown 
by MK06, therefore there is an uncertainty of the correct position 
of the extraction slit over the IFS mosaic of the order of ~ 2 arcsec. 
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solid line spectrum in the left upper panel was obtained by 
extracting from the PINGS mosaic all the individual spec- 
tra within the area enclosed by the rectangular shape and 
taking into account the overlapping regions and the covering 
fraction of the fibres due to the dithering technique applied 
to this galaxy. 

The blue line corresponds to the spectrum obtained by 
MK06 using the drift-scan technique. The dotted line was 
derived after co-adding all the spectra of the seven IFU po- 
sitions shown in Figure 11 and considering the dithering 
overlaps. This spectrum is constructed out of 6951 individ- 
ual spectra covering the whole optical B25 mag arcsec"^ 
radius of NGC 1058. In this case we can see clearly the ef- 
fect of the incomplete spatial coverage when comparing the 
integrated spectrum within the slit-aperture of MK06 and 
the whole optical radius covered by the PINGS mosaic (dot- 
ted line). Although it is beyond the scope of this paper to 
give a more quantitative comparison, we can notice from a 
purely qualitative point of view that the determination of 
some physical properties may vary when using both the slit- 
aperture and the whole mosaic integrated spectra, probably 
by as much as ~ 15-20%. 

In a similar way, the bottom panels of Figure 11 show 
the integrated spectra of the starburst galaxy NGC 3310. 
The black line shows the integrated spectrum obtained af- 
ter co-adding all the spectra within the rectangular pattern 
corresponding to the MK06 slit-aperture and considering all 
the overlaps in the mosaic due to the dithering observing 
method and the non-standard construction of the mosaic. No 
foreground stars were found within this field. The blue line 
corresponds to the MK06 drift-scan integrated spectrum. 
In this case the PINGS extraction is spatially incomplete 
when compared to the rectangular area of the MK06 aper- 
ture. Nevertheless, all the spectral features are reproduced 
quite well when we compare the PINGS spectrum with the 
MK06 data. The apparent disagreement seen in some of the 
emission lines can be explained by the spatial incomplete- 
ness of the IFS mosaic with respect to the MK06 slit, and 
the presence of a bright HII region at {Aa, AS) ~ (-10,- 
30), which is located at the edge of the extraction slit and 
presents particularly strong emission in the [O III] lines. A 
small shift of the reported position of the MK06 slit results 
in a slightly different integrated spectrum due to this strong 
H II region, causing that the peaks of the emission lines dif- 
fer in height. A stronger effect is seen in the blue spectral 
region (A < 3800 A) which is expected due to the intrinsic 
problems of the spectrophotometric calibration in this re- 
gions as explained in section 6. The integrated spectrum of 
NGC 3310 using all the fibres in the three dithered positions 
differs very little from the extracted one, and it is not shown 
for the sake of clarity. 

The spectroscopic mosaics of NGC 1058 and NGC 3310 
were not re-calibrated using the broad-band imaging tech- 
nique described in Paper II, and therefore they are not 
calibrated for the absolute zero point flux. Nevertheless, 
both the spectral shape and the spectral features of the 
PINGS extracted spectra corresponding to the slit-aperture 
of MK06 match remarkably well with these previously pub- 
lished, spectrophotometrically calibrated data, showing the 
quality of the PINGS IFS data set. 



7.2 Emission line maps 

One of the main objectives of the PINGS project is to ob- 
tain complete maps of the emission-line intensities which 
could then be analysed to describe the 2D spatially-resolved 
distribution of the physical properties of the sample. The 
PINGS spectroscopic mosaics allow us to obtain for the first 
time a complete 2D view in the optical wavelength range of 
the main emission lines used in typical abundance diagnos- 
tics methods, and important spectral features useful for the 
analysis of the underlying stellar populations. 

The ionized gas in spiral galaxies exhibits a complex 
structure associated morphologically with star-forming re- 
gions located mainly along the spiral arms. Previous at- 
tempts to perform a 2D analysis have made use of narrow- 
band and Fabry-Perot imaging at different spectral widths. 
However, in some cases the narrow-band imaging includes 
the contribution of more than one single emission line, 
such as the case of the Ha imaging, which includes the 
[Nil] AA 6548, 6584 doublet, or the [Sll] density sensitive 
doublet at AA 6717, 6731. This factor limits the utilization 
of these techniques to study just the basic parameters of 
the ionized gas, under the assumption of fixed line ratios. 
The great advantage of IFS arises from the fact that we are 
able to deblend emission lines at any discrete spatial loca- 
tion, and to ultimately produce maps of individual emission 
lines. 

In order to extract any physical information from the 
data set, we need first to identify the detected emission lines 
of the ionized gas and to decouple their emission from the 
stellar population in each individual spectrum of the IFS 
mosaic. We used population synthesis to model and sub- 
tract the stellar continuum underlying the nebular emission 
lines. This technique results in emission-line measurements 
corrected (to a first-order) for stellar absorption. The details 
of this process are described in Paper II, however we present 
briefly here the scheme followed to decouple the stellar pop- 
ulation and the emission lines in our database: i) A set of 
emission lines is identified from any strong H II region of the 
galaxy, ii) For each spectrum in the data set, the underlying 
stellar population is fitted by a linear combination of a grid 
of Single Stellar Populations (SSP), masking all the nebular 
and sky emission lines. The template models are selected in 
order to cover the widest possible range of ages and metal- 
licities. We consider the effects of dust extinction by varying 
Av from to 1 mag at A 0.2 mag. iii) We subtract the fit 
stellar population from the original spectrum to get a resid- 
ual pure emission-line spectrum, iv) Finally, we derive the 
intensities for each detected emission line. 

Individual emission-line fluxes were measured in each 
spectrum by considering spectral window regions of ~ 200 A. 
We performed a simultaneous multi-component fitting us- 
ing a single Gaussian function (for each emission line con- 
tained within each window) plus a low order polynomial 
(to describe the local continuum and to simplify the fitting 
procedure) using FIT3D (Sanchez et al. 2006). The central 
redshifted wavelengths of the emission lines were fixed and 
since the FWHM is dominated by the spectral resolution, 
the widths of all the lines were set equal to the width of 
the brightest line in this spectral region. This procedure 
decreases the number of free parameters and increases the 
accuracy of the deblending process (when required). Line 
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Figure 12. Emission line maps obtained from the IFS mosaic of NGC 1058 (without the outlying pointing). In the top-left panel, Ha 
line intensity map, top-right Ha/H/3 flux line ratio. Bottom-left, [Oil] A3727, bottom-right [OIII] A4959 -I- A5007 intensity maps. All 
maps were obtained by fitting a single Gaussian function to the emission line for each single spectrum in the mosaic. The intensity map 
is then reconstructed by interpolating the recovered flux at the location of each spectrum as described in the text. All maps arc in units 
of 10~^^ erg s~^ cm~^ arcsec"^. No correction for dust extinction was applied to any map. 



intensity fiuxes were then measured by integrating the ob- 
served intensity of each hne. The statistical uncertainty in 
the measurement of the hne flux was calculated by propa- 
gating the error associated to the multi-component fitting 
and considering the signal-to-noise of the spectral region. 

We applied this line-fitting method to NGC 1058 and 
created emission line maps by interpolating the intensities 
derived for each individual line in each individual spec- 
trum, correcting for the dithering overlapping effects. The 
interpolation was performed using E3D, adopting a natural- 
neighbour, non-linear interpolation scheme, with a final scale 
of 1" /pixel in the resulting maps. The data at the location 
of bright foreground stars in the field were masked prior to 
any interpolation, in order to decrease the effects of their 
contamination. 

As a prime example of this technique for the pur- 
pose of this paper, we calculated the following emission 
line maps: Ha, H/3, the doublet-blended [O II] A3727, and 



[O III] AA 4959, 5007. In order to prevent contamination by 
low signal-to-noise data, we masked all pixels that corre- 
spond to an integrated fiux per fibre below 10~^® erg s~^ 
cm~^. Figure 12 shows the emission line maps in units of 
10^^^ erg s~^ cm~'^ arcsec"'^. The top panels show the dis- 
tribution of the star-forming regions in the galaxy traced 
by the Hq intensity map (left) and the distribution of dust 
extinction, outlined by the Ha/H/3 fiux ratio (right). From 
these maps we note that, although many of the most intense 
emission regions are located along the inner spiral arms of 
the galaxy, some of them lie in outer regions not associated 
with any spiral structure (e.g. to the south of the galaxy, at 
ARA ~-10, ADec ~-90 arcsec). On the other hand, the 
Hq/H/3 fiux ratio shows a smooth distribution, especially 
concentrated in the inner part of the galaxy with strong 
peaks in the outlying star-forming regions. As a matter of 
fact, the very intense outlying regions of star-formation are 
responsible for the difference found in between the inte- 
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grated spectra of NGC 1058 when comparing the central slit- 
aperture of MK06 to the whole mosaic (see subsection 7.1). 
NGC 1058 is known to have a larger than average H I-to- 
optical size, and deep Hq imaging has revealed the existence 
of recent massive star formation out to and beyond two opti- 
cal radii (Ferguson et al. 1998, hereafter FGW98). However, 
the small number of H II regions analysed in previous stud- 
ies of NGC 1058 has prevented a complete explanation of 
the existence of these extreme outer regions, and of the in- 
trinsic scatter of the abundance gradient seen in this galaxy. 
Given the spatial information provided by the PINGS spec- 
troscopic mosaic, we could in principle study the properties 
of the outer HII regions with those of the inner disk and 
test the predictions of chemical evolution models and the 
behavior of the chemical abundances not only as a function 
of increasing radius, but taking into consideration the mor- 
phology of the galaxy. 

The bottom panels show the line intensity distribution 
of [Oil] A3727 (left), and [GUI] A4959-HA5007 (right). For 
both cases, we notice strong emission regions at ARA ~ - 
20 arcsec, and again several outlying emission regions that 
do not follow the spiral pattern of the galaxy. More definite 
interpretations than the qualitative ones described here can 
be achieved with a full spectroscopic analysis of the dust 
content, ionization conditions and metallicity distribution 
of the the whole IFS mosaic. 

7.3 Emission line ratios: comparison with 
literature 

Several galaxies in the PINGS sample have been spectro- 
scopically studied previously by different authors. As a con- 
sistency test of the quality of our data we performed a 
comparison of several emission line ratios as measured by 
the procedure described in the previous section with se- 
lected H II regions from the literature for which the authors 
published the emission line ratios for different species. Ta- 
ble 3 shows a comparison of the PINGS emission line in- 
tensities for three of those regions: FGW 1058E and FGW 
1058H, analysed by FGW98, and the Jumbo HII region 
in NGC 3310 observed by Pastoriza et al. 1993, (hereafter 
Pas93). 

FGW98 observed a total of 8 HII regions in NGC 1058, 
however, regions FGW 1058A to D are located in the inner 
part of the galaxy and their identification is somewhat un- 
clear (see Figure 2 from FGW98). Regions FGW1058F and 
FGW1058G fall outside the observed FOV of PINGS. On 
the other hand, FGW 1058E is a bright H II region located 
at ~ (-10, -88) arcsec in (ARA, ADec) units with respect 
to the galaxy centre (see Figure 12), while FGW 1058H is 
an outlying HII region located at ~ (140, 90) in the PINGS 
mosaic (see Figure 2). For this comparison we selected these 
two objects as both fall within the FOV observed by PINGS 
and are uniquely distinguishable from the Ha maps/images. 

In the case of FGW 1058, we used for this comparison 
a 4 arcsec circular aperture centered on the fibre with the 
strongest emission in Ha, assuming that the long-slit ob- 
servation was placed in this region as FGW98 did not give 
details of the observation of each specific object and their 
given offsets are only approximate with respect to the centre 
of the galaxy (even the extraction aperture is uncertain as 
they only quote a size range from 3 to 15 arcsec depend- 



ing on the seeing and on the size of the object in question). 
In the case of FGW 1058H, the identification was relatively 
simple as this is a small and well-defined outlying H II region. 
The emission line ratios for this region were obtained from a 
single fibre at the quoted position. FGW98 published the ob- 
served and reddening corrected emission line intensities for 
9 spectral lines, in Table 3 we compare our results with the 
observed un-corrected emission line fluxes only. Additional 
spectral lines and atomic species detected and measured in 
the PINGS data are also included in Table 3. 

Pas93 performed an optical and near-IR spectroscopic 
analysis of circumnuclear HII regions (at less than 400 pc 
from the nucleus) and the Jumbo HII region in NGC 3310. 
The identification of the circumnuclear H II regions is some- 
what difficult in the PINGS mosaic given the linear scale of 
the regions and the size of the fibres. Therefore we chose to 
analyse the bright Jumbo H II region for which emission line 
intensities were measured by Pas93. The slit position, aper- 
ture and PA are well described by these authors, however 
we did not choose to simulate an aperture in the PINGS 
mosaic as the aperture of the slit used by Pas93 (1.5 arcsec) 
is smaller than the size of a single PPAK fibre. Instead, we 
chose the fibre within this region with the strongest emis- 
sion in Ha, which corresponds to an offset (-10.6, -2.2) 
in (ARA, ADec) with respect to the galaxy centre in the 
PINGS mosaic. 

Pas93 quoted only the reddening corrected line intensi- 
ties, therefore for this comparison. Table 3 shows the ob- 
served and reddening corrected line ratios obtained from 
the PINGS spectrum extracted for this region. The PINGS 
observed line intensities listed in Table 3 were corrected 
for reddening using the Balmer decrement according to 
the reddening function of Cardelli et al. (1989), assuming 
R = Av/E{B - 1/) = 3.1. Theoretical values for the intrin- 
sic Balmer line ratios were taken from Osterbrock & Ferland 
(2006), assuming case B recombination, an electron density 
of Tie = 100 cm~^ and an electron temperature = 10* 
K. We have used only the Ha/H/3 ratio to deduce the loga- 
rithmic reddening constant c(H/?), obtaining a value of 0.32 
for the PINGS spectrum, whereas Pas93 obtained c(H/3) — 
0.29 for the same region. We did not correct the emission line 
fluxes for underlying stellar absorptions. Formal errors in the 
derived line ratios were estimated by summing in quadrature 
the uncertainty in flux calibration, the statistical error in the 
measurement of the line flux and the error in the c(H/3) term. 

Despite the ambiguity due to the exact location and 
extraction aperture of the observed spectra, from Table 3 
we can see that there is a very good agreement between the 
previously published emission line ratios and the PINGS ob- 
servations for these three H II regions. The strongest devia- 
tion is found in the [O II] A3727 as we expected, since this 
line falls in the spectral region at which the instrumental 
low sensitivity increases the flux calibration error. The com- 
parison with the NGC 1058 regions is more straightforward 
as these are the observed values without any further cor- 
rection. The comparison of the emission intensities with the 
NGC 3310 Jumbo HII region has to be made carefully since 
the values of Pas93 were corrected for underlying absorp- 
tion and using a different reddening curve. Nevertheless, the 
reddening corrected values of PINGS are in good agreement 
with the values derived by Pas93, we even detected and mea- 
sured correctly the very faint [O III] A4363 line, showing the 
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Table 3. Comparison of the PINGS emission line intensities for two HII regions in NGC 1058 and one in NGC3310 with previous 
observations and published line ratios. For NGC 1058 we compare with the regions E and H by Ferguson et al. (1998), and the Jumbo 
HII region in NGC 3310 observed by Pastoriza et al. (1993). For NGC 1058, the coordinates shown represent approximate offsets of 
the HII regions with respect to the galaxy nucleus in the PINGS mosaic in the format (ARA, ADec) with NE positive. All values for 
NGC 1058 are the observed line ratios. In the case of the NGC 3310 Jumbo HII region, the observed and redenning corrected ratios for 
PINGS are shown as Pas 93 only shows corrected line intensities (see text for details). All line ratios arc normalised to H/3. * Observed 
H/3 flux in units of 10~^® erg s~^ cm~^. 
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I(A5007)/I(A4959) 3.09 ± 0.15 2.92 ± 0.14 2.83 ± 0.20 2.85 ± 0.15 3.05 ± 0.05 3.02 ± 0.08 3.12 ± 0.03 

* F(H,3) A4861 236.2 210.5 11.6 19.5 1096 1690 



good quality of the sky subtraction despite the presence of 
the Hg A4358 sky line. The Une strength of the [O III] A4363 
line combined with the redshift of NGC 3310 makes it fall in 
the region of "detectability" as described in the simulation 
presented in subsection 6.2. 

An additional assessment of the quality of the data is 
given by the value of the I(A5007)/I(A4959) ratio which can 
be predicted from atomic theory and observed in high signal- 
to-noise astronomical spectra. According to current atomic 
computations, the theoretical value for this intensity ratio is 
2.98 (Storey & Zeippen 2000). In Table 3 we show this ratio 
for the three H II regions. The observed central values for the 
region FGW 1058E in both datasets are in good agreement 
with the theoretical value. In the case of FGW 1058H both 
values are quite similar but differ within the errors from 
the theoretical value. For the NGC 3310 Jumbo HII region, 
the central value observed by PINGS is closer to the theo- 
retical value for both uncorrected and corrected line ratios 
compared to the Pas93 value. 

The flux observed by PINGS in the H/3 hue for 
FGW 1058E is slightly higher than the one measured by 
FGW98, contrary to the case of FGW 1058H, where the flux 
in PINGS for the same line is somehow smaller, reflecting 
the unknown aperture extraction for the FGW98 long-slit 
spectrum. In the CEise of the NGC 3310 Jumbo HII region, 
the flux measured by PINGS is somewhat smaller than the 
one measured by Pas93. 



The 2D character of the PINGS data allows us to study 
the variation of the spectra within a given area that would be 
otherwise taken as a single HII region. In Figure 13 we show 
the effect of the extraction aperture on the emission line in- 
tensities for the HII region FGW 1058E. FGW98 considered 
FGW 1058E as a single H II region, however a closer look us- 
ing the dithered spectroscopic mosaic shows that this region 
is actually a complex composed of several knots and sub- 
structures with varying emission fluxes in the most promi- 
nent lines. In order to examine the difference in the emission 
line ratios in this region, we take as a central position the fl- 
bre with the strongest emission in H« within this area, with 
an integration aperture of 2.7 arcsec diameter. We then take 
concentric circular apertures of different sizes (ranging from 
4 to 12 arcsec in diameter), we integrate the spectra within 
these apertures to obtain spectra from which we measure a 
different set of emission line intensities. 

The flve top panels of Figure 13 show the variation of 
the emission line ratios obtained at different extraction aper- 
ture sizes for some relevant lines (normalised to H/3). The 
second data point corresponds to the value shown in Table 3 
(4 arcsec aperture), the horizontal line in each panel shows 
the central value and the error bar obtained by FGW98 (the 
dark/light blue colour line/band in the online version). The 
bottom panel shows the integrated flux of H/3 at each aper- 
ture. From Figure 13 we can see that the emission line ratios 
measured using different extraction apertures vary consid- 
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Figure 13. Variation of the emission line ratios as a function 
of integration aperture of the HII region FGW 1058E. The hor- 
izontal line/band in each panel shows the value/error derived by 
FGW98. All emission line intensities are normalized to H/9. The 



erably as a function of the aperture size, and that in most 
cases the dispersion of the central values is larger than the 
error of the measurements, reflecting that this is a physical 
effect. All emission line ratios tend to converge to a cer- 
tain value as the aperture size increases. Note that at the 
flux level in H/3 measured by FGW98 for this region, all the 
emission line intensity ratios, as measured by PINGS, are 
basically the same (within errors) as the values derived by 
FGW98. 

From this exercise we note that, to a first-order, the 
emission line ratios measured for a given HII region may 
significantly depend on the morphology of the region, on 
the slit (fibre) position, on the extraction aperture and on 
the signal-to-noise of the observed spectrum. All these ef- 
fects should be taken into account when deriving physical 
quantities from spectroscopic studies of H II regions. 



7.4 SN 2007gr in NGC 1058 

As another example of the potential of the PINGS spectro- 
scopic data, we present the spectrum of a particular object 
observed during the mosaicking of the galaxy NGC 1058. 
The observation of this galaxy was carried out during the 
nights of the 7th and 8th of December 2007. During this pe- 
riod we were able to observe the spectrum of the Supernova 
2007gr, previously discovered by the Lick Observatory Su- 
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Figure 14. Extracted spectrum of the SN 2007gr in the mosaic 
of NGC 1058. The flux units are IQ-^'' erg s"! cm^ A"!. Some 
nebular emission lines are seen superimposed to the SN spectrum. 



pernova Search on the I5th of August 2007 (Chornock et al. 
2007). 

SN2007gr was located at 24.8" West and 15.8" 
North of the nucleus of NGC 1058, with coordinates RA 
02'M3"" 27."98 and Dec: -f37° 20'"44."7. In the mosaic dia- 
gram of NGC 1058 shown in Figure 11, the arrow shows the 
location of SN 2007gr between two bright foreground stars. 
Figure 14 shows the optical spectrum of the supernova ex- 
tracted from the IFS mosaic, showing the typical spectrum 
of a Type Ic core collapse supernova, with the lack of hy- 
drogen, helium and silicon absorption lines. Emission lines of 
[Oil] A3727, Hq, [Nil] AA 6548, 6584 and [SII] AA 6717, 6731 
are clearly seen on top of the SN spectra, probably refiecting 
the environment of the HII region in which this supernova 
exploded. 

As described by Valenti et al. (2008), SN 2007gr showed 
an average peak luminosity but unusually narrow spectral 
lines and an almost fiat photospheric velocity profile. SN 
2007gr motivated an extensive observational campaign for 
several reasons: it was discovered at a very early stage (5 
days after the explosion), it was located in a relatively close 
distance galaxy, it was the nearest stripped-envelope carbon- 
rich SNe ever observed, and a suitable candidate for progen- 
itor search (Valenti et al. 2008). Wide-field IFS may become 
one of the main resources for SNe research groups to find 
SNe progenitors in previously observed galaxies. The pro- 
genitor supergiant stars (M > 8 Mq) that at the end of 
their lives explode as core-collapse supernovae may show a 
strong stellar spectrum that could be recovered from 2D 
spectroscopic maps with enough spatial resolution. Further- 
more, the nebular emission from the spatially adjacent spec- 
tra could provide information on the environment in which 
these SNe explosions occur. 



8 SUMMARY 

We have presented the PPAK IFS Nearby Galaxies Survey: 
PINGS, a project designed to construct 2D spectroscopic 
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mosaics of 17 nearby disk galaxies in the optical wavelength 
range. This project represents one of the first attempts to 
obtain continuous coverage spectra of the whole surface of 
a galaxy in the nearby universe. 

The final sample includes different galaxy types, includ- 
ing normal, lopsided, interacting and barred spirals with a 
good range of galactic properties and SF environments with 
multi-wavelength public data. The spectroscopic data set 
comprises more than 50000 individual spectra, covering an 
observed area of nearly 80 arcmin^. The data set will be 
supplemented with broad band and narrow band imaging 
for those objects without publicly available images in order 
to maximise the scientific and archival value of the observa- 
tions. Future plans consider to release the PINGS data set in 
order to make it freely available to the scientific community 
as a PPAK legacy project. 

The primary scientific objectives of PINGS are to obtain 
emission-line maps and moderate-resolution spectra of the 
underlying stellar population across the disks of the galax- 
ies. These spectral maps will allow us to study the spatial 
distribution of the physical properties of the ionized gas and 
the stellar components in galaxies, solving the limitations 
imposed by the small FOV and spectral coverage of previ- 
ous attempts at obtaining the 2D information of a galaxy. 
PINGS will provide the most detailed knowledge of star for- 
mation, gas chemistry and the variations across a late-type 
galaxy. This information is very relevant for interpreting the 
integrated colours and spectra of high redshift sources. The 
details provided by PINGS will contribute significantly to 
the study of the chemical abundances and the global prop- 
erties of galaxies. 

We have assessed very carefully all sources of errors and 
uncertainties found during the intrinsically complex reduc- 
tion of IPS observations. This complexity is further deep- 
ened if one considers building an IFU spectroscopic mosaic 
of a given object for which the observations were performed 
at very different stages, for some objects even spanning a 
period of years. 

To demonstrate the use of these data, we presented a 
few relevant scientific examples that have been extracted 
from the data set. We presented the integrated spectra of the 
IFS mosaics of NGC 1058 and NGC 3310 and a comparison 
with previously published data. Furthermore, we obtained 
emission line maps for NGC 1058 and we presented a quali- 
tative description of the 2D distribution of physical proper- 
ties inferred from the line intensity maps. We performed a 
comparison of the emission line ratios obtained for a number 
of selected HII regions from the literature, we showed that 
the emission line intensities obtained from the PINGS data 
agree with previously published studies, and that the line 
ratios derived from a single H II region may depend strongly 
on the position and aperture size of the slit/fibre (among 
other factors) for a given spectroscopic study. Finally, we 
showed the spectra of the SN 2007gr observed in NGC 1058 
as an example of the utility of wide-field IFS surveys to find 
potential SNe progenitors and their environmental proper- 
ties. Independent studies for individual objects, regions and 
galaxies samples will be presented in future papers. As a 
highlight of this project, the spectroscopic mosaicking of one 
of the sample objects: NGC 628, represents the largest area 
ever covered by an IFU observation. In Paper II we present 



a complete analysis of the integrated and spatially resolved 
properties of this galaxy. 
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